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^ Abstract 

o 

q ■ We present the Monte Carlo generator RacoonWW that computes cross sections 

to all processes e + e~ — > 4/ and e + e~ — > 4/7 and calculates the complete 0(a) 
electroweak radiative corrections to e + e~ — ► WW — ► 4/ in the electroweak Standard 



Model in double-pole approximation. The calculation of the tree-level processes 



D . e + e~ — > 4/ and e + e~ — > 4/7 is based on the full matrix elements for massless 

(polarized) fermions. When calculating radiative corrections to e + e~ — > WW — > 4/, 
the complete virtual doubly-resonant electroweak corrections are included, i.e. the 
factorizable and non-factorizable virtual corrections in double-pole approximation, 
and the real corrections are based on the full matrix elements for e + e~ — > 4/7. The 
matching of soft and collinear singularities between virtual and real corrections is 
done alternatively in two different ways, namely by using a subtraction method or by 
applying phase-space slicing. Higher-order initial-state photon radiation and naive 
QCD corrections are taken into account. RacoonWW also provides anomalous 
triple gauge-boson couplings for all processes e + e~ — > 4/ and anomalous quartic 
gauge-boson couplings for all processes e + e~ — ► 4/7. 
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PROGRAM SUMMARY 



Title of program: RACOON WW, version 1.3 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland 

Computer: any computer with FORTRAN 77 compiler. The program has been tested on 
DEC/ALPHA workstations and Linux PCs. 

Operating systems: UNIX, Linux 

Programming language used: FORTRAN 77 

Memory required to execute with typical data: about 4 MB 

No. of bytes in distributed program, including test data, etc.: about 2MB 

Distribution format: uuencoded compressed tar file 

Keywords: radiative corrections, W-pair production, four-fermion production, triple gauge- 
boson couplings, quartic gauge-boson couplings, Monte Carlo integration, event generation 

Nature of physical problem 

Precision calculations for W-pair production, four-fermion production, and four-fermion 
plus photon production in e + e~ annihilation. 

Method of solution 

Multi-channel Monte Carlo integration of perturbatively calculated multidifferential cross 
sections including radiative corrections to e + e" — > 4/. The code can generate weighted 
and unweighted events. For hadronization an interface to Pythia is provided. 

Restrictions on the complexity of the problem 

Virtual O{o) radiative corrections are included in double-pole approximation (DPA) only. 
Higher-order initial-state radiation is included via collinear structure functions in leading- 
logarithmic approximation. QCD corrections are included as a global multiplicative factor 
(1 + cts/n) for each hadronically decaying W boson. The external fermions are assumed 
to be massless. Below about 170 GeV the DPA is not applicable, and above 500 GeV con- 
tributions become numerically important that are not yet implemented in RacoonWW. 

Typical running time 

On a DEC/ ALPHA workstation the computing time for calculating 10 7 weighted events 
for tree-level predictions is about 2 hours. When including radiative corrections the com- 
puting time for 5 x 10 7 weighted events amounts to about 30 hours. 
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LONG WRITE-UP 



1 Introduction 

Experimental investigations of four-fermion production processes, e + e~ — > 4/(7), are 
phenomenologically very interesting for a variety of reasons. The most prominent sub- 
process proceeds via e + e~ — > WW — > 4f, which is analyzed in the LEP2 experiment 
at the per-cent level and can be studied with even higher precision at future e + e~ lin- 
ear colliders, where centre-of-mass (CM) energies go far beyond the 210 GeV reached by 
LEP2, most probably up to the TeV range. Upon reconstructing the W bosons from 
their decay products, the W-boson mass is determined at LEP2 with an accuracy of 
presently ~ 40 MeV, close to the finally aimed precision of ~ 30 MeV. Moreover, the 
triple gauge-boson couplings ZWW and 7WW are directly constrained within a few per 
cent upon analyzing angular distributions. Recently the LEP2 collaborations started to 
put bounds on specific deviations from the Standard-Model (SM) quartic gauge-boson 
couplings, mainly by investigating cross sections and photon-energy distributions in the 
processes e + e~ — > WW(7) — > 4/7. Of course, also other subchannels in four-fermion 
production, such as Z-pair or single-W production, deliver important information that is 
complementary to results obtained from W pairs, although the corresponding accuracy at 
LEP2 is significantly lower. Details about the present situation of the LEP2 data analysis 
can be found in Ref. [1] and about future prospects at linear colliders in Ref. [2,3]. 

To account for the high experimental accuracy on the theoretical side is a great chal- 
lenge: the W bosons have to be treated as resonances in the full four-fermion processes 
e + e~ — > 4/(7), and radiative corrections need to be included. From the variety of the 
above-mentioned physical issues it is also clear that semianalytical predictions are not 
sufficient; the data analyses require realistic studies with event generators that meet 
these theoretical requirements. In this article we describe the structure and the use of 
the Monte Carlo (MC) event generator RacoonWW, which is designed for this task. 
Numerical results of RacoonWW as well as the underlying calculation of the cross sec- 
tions and corrections have already been published in Refs. [4,5,6,7,8]. A general review of 
the state-of-the-art for four-fermion production can be found in Refs. [9] and [10]. The 
latter reference, in particular, contains detailed tuned comparisons of RacoonWW re- 
sults to related calculations, such as the semi- analytical calculation of Ref. [11] and the 
MC generator YFSWW3 [12] for e + e~ — > WW — > 4/ including 0(a) corrections in 
double-pole approximation (DPA), and the MC generators WRAP [13] and PHEGAS 
[14] for tree-level predictions of e + e~ — > 4/7. 

RacoonWW operates in three different modes: the tree-level mode, the improved Born 
approximation (IB A) mode, or the radiative- correction (RC) mode. The tree-level mode 
provides predictions for e + e~ — > 4/ or e + e~ — > 4/7 for all possible final states with 
massless fermions in lowest order, based on the full set of Feynman diagrams or on selected 
subsets [4] . Gauge-boson decay widths are implemented as constant or running widths, or 
within the complex-mass scheme. Since the latter fully respects gauge invariance, gauge- 
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invariance violations in the other schemes can be controlled. The initial-state radiation 
(ISR) in leading-logarithmic (LL) approximation via the structure-function method and 
the Coulomb singularity for slowly moving W bosons can be optionally included [7,8]. 
Anomalous triple gauge-boson couplings are supported for all processes e + e~ — > 4/ [15], 
and for the process class e + e~ — » 4/7 anomalous quartic gauge-boson couplings are also 
available [8]. 

The IBA mode calculates the leading universal corrections to the e + e~ — > WW — > 4/ 
processes [7] or to the process class e + e~ — » 4/7 [8]. The IBA for e + e~ — ► WW — > 4/, 
which is similar to the approach followed in the semi- analytical program GENTLE [16], 
is particularly interesting near the W-pair production threshold where it is still applicable 
in contrast to the full one-loop calculation in DPA. 

The RC mode, which is designed for e + e~ — > 4/(7), uses the 4/ lowest-order contribution 
of the tree-level mode and further includes the full electroweak 0(a) corrections in the 
DPA as defined in Ref. [6]. More precisely, only the virtual corrections are treated in the 
DPA leading to factorizable and non-factorizable corrections [17]. The former comprise 
the corrections to the production and decay subprocesses, while the latter are due to 
soft-photon exchange between these subprocesses. The real 0{a) corrections are based 
on the smallest gauge-invariant subset of the full e + e~ — > 4/7 matrix element that in- 
cludes W-pair production. They are evaluated analogously to the tree-level mode for 4/7 
production. Since the full calculation is consistently based on massless external fermions, 
photons that are collinear to any external charged fermion have to be recombined with 
this fermion, as described in Ref. [6] in more detail. RacoonWW provides such recombi- 
nation algorithms, which can, however, be replaced by user-defined alternatives. Finally, 
higher-order ISR at the LL level is included up to 0(a 3 ) using the structure-function ap- 
proach [18]. In this context, it should be emphasized that the combination of the virtual 
corrections in DPA with the real corrections based on full matrix elements is non-trivial. 
The singularity structures of both parts are closely related, and a mismatch has to be 
carefully avoided. RacoonWW provides two versions for the redistribution of the soft 
and collinear singularities, resulting in two different branches, called slicing and subtrac- 
tion branches. The former employs the well-known phase-space slicing, the latter the 
subtraction procedures of Refs. [19,20]. Since even the phase-space generation in the two 
branches is different, a comparison of results obtained with the different branches is a 
highly non-trivial consistency check. 

Finally, RacoonWW offers various further options, such as the choice between differ- 
ent input-parameter schemes, the inclusion of naive QCD corrections, a switch for the 
imaginary parts of loop diagrams, and various settings to estimate the intrinsic uncer- 
tainty of the DPA. RacoonWW is able to generate unweighted events by a hit-and-miss 
algorithm. In the RC mode unweighted events can only be generated when using the 
slicing branch. In order to account for the complex peaking behaviour of the matrix ele- 
ments RacoonWW uses a multi-channel integration technique [21] and adaptive weight 
optimization [22]. 

RacoonWW has been tested up to a CM energy of y/s — 500 GeV in the RC mode. 



4 



The theoretical uncertainty of the total cross section has been estimated to be 0.4% 
for energies between 200 GeV and 500 GeV, 0.5% at 180 GeV, and 0.7% for 170 GeV 
[6,10]. The estimates result from investigations of intrinsic ambiguities in the definition 
of the DPA and from a detailed comparison with results obtained from the MC generator 
YFSWW3 [12]. In the RC mode of Racoon WW results have also been obtained for a 
CM energy of 1 TeV, but no detailed study of the theoretical uncertainties has been per- 
formed. Around the W-pair production threshold at 160 GeV the DPA is not applicable, 
and above 500 GeV higher-order radiative corrections that are not yet implemented in 
Racoon WW, e.g. large logarithmic electroweak corrections beyond 0(a), may become 
numerically sizeable. In the tree-level mode RacoonWW is tested up to CM energies 
of 10 TeV [4], and IBA results have been obtained for CM energies up to 1 TeV [7]. The 
tree-level and IBA modes are also applicable in the W-pair production threshold region. 

The paper is organized as follows: in Section 2, we sketch the general structure of the 
program. Section 3 is devoted to the description of important details of the code. In 
Section 4, we describe in detail the usage of RacoonWW including all important options. 
Section 5 summarizes the paper, and the appendix contains input and output of two 
sample runs. 



2 Program structure 

RacoonWW provides predictions for the cross section to the class of processes e + e~ — > 
4/(7) as given by 

/da = M^ Jd*™ \M( P+ ,p-,h,...,k n+2 )\ 2 , (1) 

where p± and ki are the momenta of the incoming and outgoing particles, respectively, 
and n = 4, 5 is the number of outgoing particles. The total energy of the scattering process 
in the CM frame, where p + + p = 0, is denoted by y/s = p° + +p°_. The n-particle phase 
space is defined by 

™ +2 H 3 k. / ™+ 2 \ 

'"i V 1 " v i 

Results for the matrix element squared \Ai\ 2 in terms of helicity amplitudes for the tree- 
level processes e + e" — > 4/(7) with massless fermions are given in Ref. [4], and the cal- 
culation of radiative corrections is described in Ref. [6]. The phase-space integration is 
performed with the help of a MC technique, which allows to calculate a variety of observ- 
ables simultaneously and to easily implement separation cuts and photon recombination 
schemes in order to account for the experimental situation. The matrix element squared 
exhibits a complicated peaking behaviour in various regions of the integration domain. 
In order to obtain a numerically stable result and to reduce the MC integration error, 



(2) 
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Racoon WW is based on a multi-channel MC method [21] with adaptive weight opti- 
mization [22]. 

The general structure of the program is determined by the generation of the n-particle 
phase space, and the evaluation of the matrix element squared, \M\ 2 . The program 
run starts with some initializations and ends with the output of the integration result in 
form of an output file for the total cross section and a number of data files for various 
distributions. 

The initialization of the matrix-element evaluation consists to a large extent of specifying 
the input parameters and the different options provided by RacoonWW. The initializa- 
tion of the phase-space generation in the multi-channel MC integration approach essen- 
tially amounts to the identification of kinematical channels for the given final state. While 
the user is free to choose the parameters and flags for options in the input file (and to 
a lesser extent also in racoonww.f and public. f), the initialization of the multi-channel 
phase-space generators is done automatically by RacoonWW. 

RacoonWW consists of 11 fortran files and a makefile. The usage of the program is 
described in the file README. RACOONWW also includes 5 default input files with corre- 
sponding output and data files. They represent the recommended choices for input param- 
eters and flags of the different modes provided by RacoonWW. The general structure of 
the program is reflected by a series of subroutine calls in the main program in racoonww . f : 

CALL READ INPUT : the input file is read. 

CALL INITIALIZE : parameters and flags are set, and the MC integration is ini- 

tialized (specification of channels for multi-channel integra- 
tion). 

CALL INTEGRATION : the MC integration and the generation of NEVENTSW weighted 

events are performed, the hit-and-miss algorithm for gener- 
ation of unweighted events is initialized. 

CALL EVENTGENERATION : in the main program it is only called if NEVENTSUNW > 0; 

then NEVENTSUNW unweighted events are generated by using 
a hit-and-miss algorithm (see Section 4.4). This subroutine is 
also called in the subroutine INTEGRATION where it generates 
weighted events. 

CALL PY4FRM : the interface to PYTHIA [23] to include parton showering and 

hadronization is called. The subroutine PY4FRM is provided 
by Pythia (see also Section 4.4). 

CALL WRITEOUTPUT : the results are written to the output file and data files. 

Besides the choices made for the input parameters and flags in the input file, the user can 
change the default values for the remaining parameters and flags of RacoonWW in the 
subroutine INITIALIZE in racoonww.f and in the subroutine PARAMETER in public. f. 
The specification of the input in Racoon WW is described in detail in Section 4.1. After 
the initialization, the MC integration is performed in 
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SUBROUTINE INTEGRATION (AVPREL , SIPREL) 

AVPREL : total cross section calculated with NEVENTSW weighted events, 
SIPREL : corresponding statistical error 

by calling the MC kernels in 

SUBROUTINE EVENTGENERATI ON (ATOTSQR , A1SQR , A2SQR , UNWEIGHTING) 

ATOTSQR : total matrix element squared for the event (input for PY4FRM) , 

A1SQR : matrix element squared for the configuration with fermions fs, f± and 

/s, f§ as the two colour singlets (input for PY4FRM), 

A2SQR : matrix element squared for the configuration with fermions fa, f$ and 

/ 5 , /4 as the two colour singlets (input for PY4FRM), 

UNWEIGHTING : flag for generating weighted (UNWEIGHTING = 0) or unweighted 
(UNWEIGHTING = 1) events; this flag is automatically controlled by 
NEVENTSUNW. 

RacoonWW provides an unweighting procedure based on a hit-and-miss algorithm. It 
is activated when the user sets NEVENTSUNW > in the main program, where NEVENTSUNW 
is the number of unweighted events. When generating unweighted events, parton show- 
ering and hadronization for coloured final states can be taken into account by activat- 
ing the interface to PYTHIA [23] in the main program. The information about the (un- 
weighted) event is stored in the HEP standard format [24] in the common block HEPEVT 
(see Section 4.4). The variables ATOTSQR, A1SQR, and A2SQR are provided for the genera- 
tion of unweighted events. A detailed discussion of the unweighting procedure is given in 
Section 4.4. 

By default RACOON WW generates only weighted events (NEVENTSUNW = 0). The number 
of weighted events, NEVENTSW, is set in the input file. The events are represented by 
the four- momenta of the final-state particles and their corresponding weights. A weight, 
w(ro,r- ? ), assigned to the MC point (rj), r J ') of the jth event, is composed of the total 
density of the event, g to t, (cf. Section 3.1.2) and the matrix element squared as 

w{ri^) = with f{h) = \M{h)\\ (3) 



where Tq, r J is a set of pseudo-random numbers. As usual, the MC estimate I of the total 
cross section is obtained as (N = NEVENTSW) 



1 N 



N 



J2w{r J ,r j ). (4) 



In the RC mode of RacoonWW more than one weight per (weighted) event is generated 
corresponding to the kinematics of the subcontributions, i.e. whether the weight is a 
4/ or 4/7 event with or without ISR, or corresponding to the DPA part of the virtual 
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corrections, or to one of the subtraction terms. The weights for each weighted event (see 
Tables 5-7) together with the corresponding four-momenta are provided in the common 
block EVENT in racoonww.f: 



COMMON/EVENT/WEIGHT , PBEAM , P , WEIGHTBORN , PBORN , NSUBEVENT 
WEIGHT(1:29) : all weights generated for one weighted event, 
PBEAM (I, J) : four-momenta of incoming particles (before ISR), 
PBEAM(1, 0) + PBEAM(2, 0) = y/s, 
1 — 1,2 : incoming positron, electron, 
J = 0, . . . , 3 : four- vector components, 
P(I,J,K) : four-momenta (after possible ISR), 
P(1,0,K) + P(2,0,K) < yfe, 
I — 1, ... ,7 : numerates the external particles 

(1 = incoming positron, 2 = incoming electron, 
3,. . . ,6 = outgoing fermions, 7 = photon), 
J = 0, . . . , 3 : four-vector components, 
K = 1, ... ,29 : numerates the weights, 
weight for tree-level four-fermion process, 
momenta corresponding to WEIGHTBORN, 
I = 1, . . . , 6 : numerates the external particles, 
J = 0, . . . , 3 : four- vector components, 
number of weights corresponding to one weighted event, 
SMC =1,3:8 weights per event (slicing branch), 
SMC = 2 : 29 weights per event (subtraction branch). 



WEIGHTBORN 
PBORN (I, J) 



NSUBEVENT 



These two key components of the event generation, phase-space generation and matrix- 
element evaluation, are described in detail in Sections 3.1 and 3.2, respectively. 

For a realistic simulation of e + e~ — > 4/(7) processes separation cuts and a photon re- 
combination procedure have to be applied to take into account the detector response. 
RacoonWW provides several experiment-inspired recombination procedures and sets of 
separation cuts, which can be specified in the input file. They are defined in public.f in 



SUBROUTINE RECOMBINATION (PBEAM , PIN , POUT , N0CUT) 

PBEAM (1 :2,0:3) : four-momenta of incoming particles (before ISR), 



PIN(1:7,0:3) 

P0UT(1:7,0:3) 

N0CUT 



four-momenta of external particles (input), 

four-momenta after eventual recombination (output), 

indicates whether the event is rejected (N0CUT = 0) or has passed 
all cuts (N0CUT = 1) 



and 
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SUBROUTINE CUT (PBEAM , P , NPARTICLE , NOCUT) 

PBEAM(1 :2,0:3) : four-momenta of incoming particles (before ISR), 
P(l:7,0:3) : four-momenta of external particles (after possible ISR), 
NPARTICLE : number of particles in the final state, 

PBEAM, NOCUT : see subroutine RECOMBINATION. 

Finally, the data files for the histograms are filled in subroutine SETTING in public. f, 
and the results for the total and differential cross sections are written to the output and 
data files in the subroutine WRITEOUTPUT. 

In the following we describe in detail the most important features of the program and its 
usage. 

3 Details of the program 

RACOONWW is a merger of two MC programs, which use the same routines for the 
evaluation of the matrix elements but different, completely independent routines for the 
phase-space generation. Moreover, they differ in the treatment of the soft and collinear 
singularities arising in the calculation of the radiative corrections to e + e~ — * WW — > 4f. 
One employs the subtraction method of Refs. [19,20] and the other applies phase-space 
slicing (for a review see e.g. Ref. [25]). The two branches of Racoon WW are named 
accordingly, the subtraction and slicing branches, and the corresponding MC integration 
kernels are 

FUNCTION KERN(XV,WGT) 

XV (1 : 15) : pseudo-random numbers for the generation of one weighted event, 
WGT : not used in this version of RACOONWW 

for the slicing branch in kern . f and 

SUBROUTINE SUBTRACTION (RANDOM) 

RANDOM (1 : 15) : pseudo-random numbers for the generation of one weighted event 

for the subtraction branch in subtraction. f, respectively. The pseudo-random number 
generator RANS can be found at the end of racoonww.f. The user chooses between these 
two branches by setting the flag SMC in the input file, SMC = 1 or 3 for the slicing branch and 
SMC = 2 for the subtraction branch. The value SMC = 3 has to be used when unweighted 
events are generated in the RC mode. The initialization of the event generation and the 
setting of branch-specific options for the slicing and subtraction branches of RacoonWW 
is done in 

SUBROUTINE INITSLICING (ENERGY , CONV , SCALE , NEVENTSW , NOUT , I0UT) 
ENERGY : CM energy in GeV (ENERGY = y/s), 
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CONV 



: conversion factor from GeV 2 to fbarn, 



SCALE : QED splitting scale in the structure functions (see Ref. [6] and Section 3.2), 
NEVENTSW : number of weighted events, 



in slicing. f and 

SUBROUTINE INITSUBTRACTION (ENERGY , CONV , SCALE , NOUT) 
ENERGY, CONV, SCALE ,N0UT : see subroutine INITSLICING 

in subtraction. f. In particular, the specification of the kinematical channels for the 
multi-channel MC integration for the chosen four-fermion final state is done automatically 
in these routines. 

After the initialization, the MC kernels are called. The phase-space cuts are applied first, 
and if the event is not discarded, the total phase-space weight, l/ptot, and the matrix 
element squared, \M.\ 2 , of Eq. (3) are calculated. For the matrix-element evaluation the 
user can choose between the tree-level, IBA, and RC modes of RacoonWW by setting 
the flag SRC in the input file. The implementation of the three modes is described in detail 
in Section 3.2. 



3.1 The phase-space generation and multi- channel MC integration 

The numerical integration of Eq. (1) is rather complicated owing to the rich peaking 
structure of the integrand, since the amplitudes to e + e~ — * 4/ and, in particular, to 
e + e~ — > 4/7 involve a huge number of propagators that become resonant or are enhanced 
in various phase-space regions. To obtain reliable numerical results, events have to be 
sampled more frequently in the integration domains where the integrand is large. This 
redistribution of events is called importance sampling. In practice, this means that the 
mapping of the pseudo-random numbers < rj < 1, into the space of final-state 
momenta has to be chosen such that the corresponding Jacobian 1/ g(r) cancels the prop- 
agators of the differential cross section at least partially. 



3.1.1 Propagator mappings 

In order to smoothen a single propagator with momentum transfer q in the square of a 
matrix element, 



NOUT 



I OUT 



: number assigned to the output file, 

: number assigned to the file optimization, info 



|M| 2 oc 



1 



(5) 



(g 2 -M 2 ) 2 + M 2 r 2 ' 
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the phase space is parametrized in such a way that the virtuality q 2 of the propagator 
is chosen as an integration variable. The integral over q 2 is transformed into an integral 
over the random number r: 



<7 2 
imax 



/ d,2 = / dr W?W)' (6) 



The dependence of q 2 on the pseudo-random number r has to be chosen such that the 
Jacobian l/g pTop smoothens the propagator (5). In RacoonWW two types of mappings 
are used according to the mass M of the propagating particle. Since all occurring massive 
particles possess a finite decay width T, a suitable mapping is of Breit-Wigner type 

(M^o, rvo): 



q 2 {r) =Mrtan[y 1 + (y 2 - Vl )r] + M 2 , 
MY 

9*™M 2 ) = {y2 _ yi)[{q2 _ M2)2 + M2 p 2] (7) 



with 



(q 2 - -M 2 \ 

j/1,2 = arctan I '-^ I . (8) 



For massless particles the following mappings are appropriate (M — F — 0): 



q 2 (r) 



# P rop(g 2 ) = 



1 - v 



1-1/ 



i 



2 Nl-i/ 
max/ 



\HmmJ J 



(9) 



? 2 W = (?Lx) r (?Ln) 1 r > Vop(g 2 ) = 2] -7 2 7-1 v - (10) 

The mapping of Eq. (9) is only valid for v ^ 1, while Eq. (10) is a substitute of Eq. (9) 
for v — 1. The naive expectation z/ = 2 from the squared matrix element is not necessarily 
the best choice because the propagator poles are partially cancelled in the collinear limit. 
It turns out that a proper value is v < 1. 

In this way all propagators of a given Feynman diagram can be smoothened by choosing 
an appropriate set of integration variables and mappings as given in Eqs. (7)-(10). The 
integration variables that do not correspond to a specific propagator mapping are sampled 
uniformly, i.e. 

Q 2 (r)=q 2 milx r + q 2 m . m (l-r), 0(r) = 2vrr, cos£(r) = 2r-l (11) 
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for invariant masses q 2 , azimuthal angles 0, and polar angles 9, respectively. 



3.1.2 Multi-channel approach 

Obviously, importance sampling of all propagators appearing in the amplitude is not 
possible by a single phase-space mapping. The multi- channel approach [21] provides a 
solution to this problem. For each peaking structure a suitable set of invariants is chosen, 
and accordingly a mapping of random numbers r into the final-state momenta fcj, so that 
the resulting Jacobian of this mapping 1/gk cancels this particular peaking behaviour of 
f(ki). The phase-space integral of Eq. (1) reads 

r ) Nch r 1 

/ = / dr £ 9(r - fa^)9(fa - r ) / d 3 "~ 4 r - (12) 

J { k =i { 9tot{ki{r)) 



where the fa define a partition of unity, i.e. fa = J2)=i a j, fa — 0, /3jv ch = J2f=i a j — 1, 
and iV c h denotes the number of all channels. In Eq. (12) an additional random number 
r is introduced in order to select a channel k randomly with probability ak > 0. The 
numbers are called a priori weights. All local densities gu are combined into the total 
density g tot , 

N ch 

9tot = Y; a k9k, (13) 

k=l 



which is expected to smoothen the integrand over the whole phase space. In this way, it is 
possible to simultaneously include different phase-space mappings which are suitable for 
different parts of the integrand. To be specific, for each (tree-level) Feynman diagram one 
channel exists, so that all propagators (except massive t-channel propagators in the slicing 
branch) of the squared diagram are smoothened by the corresponding local density g^. 
The local densities are obtained as the product of the corresponding Jacobians, as given 
in Eqs. (7)-(10). 

The a priori weights are adapted in the early phase of the MC run several times to 
optimize the convergence behaviour of the numerical integration. This adaptive weight 
optimization [22] increases the a priori weights for such channels that give large con- 
tributions to the total cross section. After a certain number of generated events a new set 
of a priori weights a£ cw is calculated according to 



N 



a k oc ak 



-Y 



g k (kj(^))[w(r 3 ,Ti)f 
0tat(fci(r J ')) 



N ci 



J2 a 

k=l 



new i 

k — l - 



(14) 
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The error of the MC integration can be estimated as 

ST = \j W{a) N ~ P , W{a) = 1 J2[w(ri, r>)? (15) 

with the MC estimate I of Eq. (4). 

The phase-space integration is implemented in RacoonWW as follows: 

• First the phase-space generators are initialized by identifying all channels needed to 
account for the peaking behaviour of the process under consideration. For the processes 
e + e~ — > 4/, there are between 6 and 128 different channels, for e + e~ — > 4/7 between 
14 and 928 channels. Each channel smoothens a particular combination of propagators 
that results from a characteristic Feynman diagram. The identification of the channels 
in the slicing and subtraction branches of Racoon WW is described in Section 3.1.4 
and Section 3.1.5, respectively. 

• The channel that is used to generate the event is picked randomly with probability a^. 
Using the randomly chosen channel, the momentum configuration ki is determined from 
the corresponding phase-space generator. Finally, using the momentum configuration 
ki the local densities of all channels g^ are calculated. From those densities the total 
density g tot of the event is determined. Apart from the evaluation of the matrix elements, 
the calculation of the local densities is the most time-consuming part of the numerical 
integration. 

• In order to reduce the MC error further, the a priori weights a fe are recalculated several 
times in the MC integration according to Eq. (14). 



3.1.3 The phase-space generators 

The actual calculation of the event kinematics is decomposed into three steps: the calcu- 
lation of time-like invariants s,, of 2 — * 2 scattering processes, and of 1 — > 2 decays. The 
phase-space integration of Eq. (1) reads 

n-2 Si 7 ax r n-l-r 

/d$w=n / ^.n d$r 2 n / d< ^ 2 > 

J i=l s . . j=l J k=l 

where the phase spaces of the 2^2 scattering processes and the 1 — > 2 decays are 
denoted by &*r* 2 and $^ 2 , respectively. The n — 2 invariants Sj result from the phase- 
space factorization into scattering processes and decays. The number r of 2 — > 2 scattering 
processes can range from 1 to n — 1. If there is only one 2-^2 scattering process (r = 1), 
this can involve either a t-channel or an s-channel propagator. If there are more scattering 
processes (r > 2), these involve only t-channel propagators. If only one t-channel 2^2 
scattering process is present, i.e. for the first two topologies in the first row of Figure 1 
and the first three topologies in the second row of Figure 1, all the correspond to the 
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S2 



S3 
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Fig. 1. Generic diagram of a specific channel 

virtualities of propagators and are mapped as described in Section 3.1.1. The angles and 
invariants appearing in the parametrizations of the phase spaces are sampled as described 
below. If several t-channel propagators are present, i.e. for r > 2, some of the Sj do not 
correspond to virtualities of propagators in the diagram; for such variables no mappings 
are introduced, they are generated uniformly. Topologies with no t-channel propagator, 
are obtained from those with only one t-channel propagator by replacing the t-channel 
scattering process with an s-channel scattering process. The corresponding invariant is 
fixed to the CM energy squared s, and the phase space is parametrized as for the 1 — ► 2 
decays in Eq. (20) or equivalently for the 2 — > 2 scattering processes in Eq. (18) without 
propagator mapping (u = M = T = 0). Diagrams with four-particle vertices result from 
the contraction of propagators in the diagrams of Figure 1. Consequently, they also fit into 
the decomposition (16), and the invariants corresponding to the contracted propagators 
are sampled uniformly. More details on phase-space parametrizations can be found in 
Refs. [20,26]. 

For each class of diagrams with the same propagator structure we adopt an own channel for 
the numerical integration. No special channels are provided for interference contributions. 



We illustrate the general strategy by considering the specific topology shown in Figure 1, 
where a particular choice for the incoming and outgoing momenta p± and hi is made. For 
this example the phase-space parametrization Eq. (16) reads 



/d$ (n) | F . gi = n / ds t Jd^ 2 ( P+ ,p^h,k U57 ) Jd^ 2 (p_, P+ -h;k 3 ,h 57 ) 

i=l c ■ 

i , mm 

x J d^ 2 ( P+ - k 6 , v- - h; h, h 7 ) J d^ 2 (k 4} k 7 ), (17) 



where sums of outgoing momenta are abbreviated by fcy = hi + kj, = fcj + kj + fc* 
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etc. Note that only the time-like invariant s\ = k\ 7 corresponds to a virtuality of a 
propagator in the diagram, while s 2 = k\ 57 and s 3 = A:f 457 correspond to invariant masses 
of fictitious final-state particles within the first two 2 — > 2 scattering processes, symbolized 
by the arches in Figure 1. More examples for the phase-space generation can be found in 
Refs. [4,20]. 

In detail we proceed as follows. 



Calculation of time-like invariants 

First of all, the time-like invariants Sj are determined. The invariants corresponding to 
s-channel propagators are calculated from Eqs. (7)-(10). All other time-like invariants are 
sampled uniformly, i.e. calculated from Eq. (11). The determination of the invariants Sj 
is ordered in such a way that the Sj nearest to final-state particles are calculated first, 
followed by the next-to-nearest s-channel propagators, and so on. The lower and upper 
limits s^min/max! m general, are functions of the already determined invariants. To enhance 
the efficiency of the numerical integration in the presence of separation cuts, in some cases 
these cuts are taken into account in the evaluation of s^mm/max- 



2^2 scattering processes 

The phase spaces for the 2^2 scattering processes of Eq. (16) are parametrized as 

, 27T ^max 

/ d^ 2 ( qi ,q 2 ; Pl ,p 2 ) = I d0* / dt, (18) 

J 8J( gi q 2 y - q\q\ { ^ 



where q± :2 are the incoming and pi j2 are the outgoing momenta, and <p* is the azimuthal 
angle defined by p\ and q\ in the CM frame of the subprocess where qi + <\ 2 = 0. The 
argument t = (pi—qi) 2 of the t-channel propagator is calculated according to Eqs. (7)-(10) 
for t-channel scattering processes (apart from massive t-channel propagators in the slicing 
branch that are sampled uniformly) and according to Eq. (11) for s-channel scattering 
processes. The azimuthal angle 0* is sampled uniformly. Since the corresponding polar 
angle cos 9* of this scattering process depends on the invariant t linearly, 



_ q\p\ + q\p\ - 2(q 1 q 2 )(p 1 p 2 ) + 2^J( qi q 2 ) 2 - q\ q\ yWs) 2 - vlvl cos 9* 

t ~ ? — ; — ^2 ' 



the calculation of the momenta and the Lorentz transformation into the laboratory (LAB) 
frame is simple and, thus, skipped here. Explicit formulas can, for instance, be found in 
Refs. [20,26]. 

Cuts on angles between outgoing particles and the beams effectively reduce the integration 
range of t. For the 2^2 scattering processes that are attached to the original incoming 
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Fig. 2. The generic phase-space generators symbolized by the topologies of the corresponding 
diagrams. The arcs indicate the intermediate 2 — > 2 scattering processes used in the phase-space 
parametrizations. Topologies with an s-channel 2 — > 2 scattering process are obtained from the 
ones with only one t-channel scattering process by contracting the t-channel propagator. 

particles these cuts are in some cases included in the determination of t m in/max to enhance 
the efficiency. 



1 — > 2 particle decays 

The phase-space integration for the 1 — > 2 particle decays in Eq. (16) reads 



J d^ 2 ( Pl ,p 2 ) 



(pm) 2 - p \p\ 

4(pi + p-i, 



2„2 2n 



d(j)* / dcos0*, 



(20) 



-i 



where 0* and 6* are the azimuthal and polar angles in the rest frame of the decaying 
particle, respectively. The variables <fi* and cos 9* are sampled uniformly. The momenta of 
the outgoing particles are pi and p^- As in the former step, the calculation of the momenta 
and the Lorentz transformation into the LAB frame is straightforward. 

All topologies corresponding to the generic phase-space generators for the processes 
e + e~ — > 4/ and e + e~ — > 4/7 are shown in Figure 2. 



3.1.4 Phase-space generation in the slicing branch 

In the slicing branch, the 2 — > 2 scattering processes and the 1 — > 2 decays are combined 
in the function 

d$ 1->2 (pi,p 2 ) for MAPI, MAPIG = 0, 

# 2 (sj,MAPl/MAPIG) = { \ (21) 

d^ 2 (q 1 ,q 2 ;p 1 ,p 2 ) for MAPI, MAPIG = 1, 2, 3, 
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MAPI(6) = 1: 



MAPI(6) = 2: 




81 



S2 




MAPIG(6) = 1: 



MAPIG(6) = 2: 



MAPIG(6) = 3: 




Fig. 3. The basic phase-space decompositions from which all channels for the processes 
e + e~ — ► 4/ and e + e~ — > 4/7 can be constructed, as used in the slicing branch 

where Sj = (pi+P2) 2 - Using the function R 2 , the phase-space integration of Eq. (16) reads 



/ 



d$W= / dsi / ds 2j R 2 (si,MAPl(5)) R 2 {s 2 , MAPI (4)) R 2 (s, MAPI (3)), 

d$ (5) = J dsi J ds 2 J ds 3j R 2 (si,MAPIG(5)) 

x R 2 (s 2 , MAPIG(4)) R 2 (s 3 , MAPIG(3)) R 2 (s, MAPIG(8)). 



(22) 



(23) 



The propagator mappings for the functions R 2 and for the variables s±, s 2 , S3 are specified 
by the values of the variables MAPI and MAPIG for the four- and five-particle phase spaces, 
respectively. In this way all channels for the multi-channel integration of the processes 
e + e~ — ► 4/ and e + e~ — > 4/7 can be constructed from the two and three basic phase-space 
decompositions, respectively, shown in Figure 3. Note that several topologies of Figure 2 
fit into single structural diagrams of Figure 3. The definitions of the variables Sj and their 
phase-space boundaries are given in Table 1. 

The phase-space generators of individual channels are built generically from these basic 
phase-space decompositions in kern.f by calling subroutines PHSPGEN_4 or PHSPGEN_5: 

SUBROUTINE PHSPGEN_4 ( IDENS , II , MAPI , VZI , NU , S , X , WPS , ICUT) 

SUBROUTINE PHSPGEN_5 ( IDENS , I IG , MAPIG , VZIG , NU , S , X , WPS , ICUT) 

IDENS : determines whether the four-momenta ki (IDENS = 0) or the density 

is calculated (IDENS = 1), 

II,IIG : II (1 : 4) ,IIG(1 : 5) specify the permutations of the final-state particles, 
II, IIG = 1, . . . , 4 denote the final-state fermions in the order of the input 
file, IIG = 5 denotes the photon, 
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Table 1 

The definition of the invariants Sj in Eqs. (22) and (23) and their phase-space boundaries for 
the different topologies in the slicing branch 

MAPI, MAPIG : MAPI (6)= 1,2 and MAPIG(6)= 1,2,3 numerate the basic topologies, 



MAPI (1:2), MAPIG (1:2), and MAPIG (7) specify the s-channel propagator 
mappings of the variables s±, s 2 , S3 of Eqs. (22) and (23): 

MAPI, MAPIG = : no mapping, 

MAPI, MAPIG = 1 : mapping for a massless propagator Eq. (9) or Eq. (10), 

MAPI, MAPIG = 2 : mapping for a W-boson propagator Eq. (7), 

MAPI, MAPIG = 3 : mapping for a Z-boson propagator Eq. (7), 

MAPI (3:5), MAPIG (3 : 5) , and MAPIG (8) specify the t-channel propagator 
mappings in the functions R 2 of Eqs. (22) and (23): 

MAPI, MAPIG = : no mapping (s-channel 2^2 scattering), 

MAPI, MAPIG = 1 : no mapping (massive particle in t-channel propagator), 

MAPI, MAPIG = 2 : mapping for a massless t-channel propagator 

using Eq. (9) or Eq. (10), 
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VZI, VZIG : VZI(1:2) ,VZIG(1:3)= ±1, the final-state momenta of the corresponding 





2^2 scattering process are interchanged if VZI, VZIG = — 1, 


NU 


parameter v of Eqs. (9), 




i^ivi eneigy scjuaieQ, 


X 


3n — 4 random numbers, 


WPS 


local density of the channel, 


I CUT 


the event is rejected if ICUT = 1. If separation cuts are directly taken into 
account in the phase-space generation (see e.g. Eq. (A. 4) of Ref. [4]), it 
can happen that s min > s max or i min > t max . Then ICUT is set to 1. 



The channels for 4/ and 4/7 production are obtained from the basic topologies by speci- 
fying the parameters II, MAPI, VZI and IIG,MAPIG,VZIG, respectively. The parameters 
VZI, VZIG distinguish between t-channel processes qi + q 2 — > P\ + P2 with invariants 
t — q±—pi or t = q± — pi- For instance, the topology of Figure 1 where the final-state par- 
ticle with momentum ke is emitted from the electron instead of the positron is obtained 
by changing the sign of VZIG(l). The specifications of the relevant parameters for each 
channel for a given final state is done in the subroutine ADDMAP in slicing. f. For the 
example of the multi-peripheral diagram given in Ref. [4] , the multi-channel parameters 
are MAPIG(1:8)= {1,3,2,0,0,3,0,2}, IIG(1:5)= {3, 5, 4, 2, 1}, and VZIG(1:2)= {1,-1} 
(VZIG (3) is not used for MAPIG(4) = topologies). 

The propagator mappings specified by the parameters MAPI and MAPIG are performed 
according to Eqs. (7)-(10) in kern.f by the subroutines SMAP and TMAP, respectively. The 
corresponding Jacobians (WPSI) are calculated in the subroutines SDENS and TDENS. The 
Lorentz transformations used in the calculation of the outgoing momenta consist of a 
combination of boosts along the z direction and rotations described by the subroutines 
B00ST_Z and ROTATION in kern.f. 

3.1.5 Phase-space generation in the subtraction branch 

In the subtraction branch, the generation of the momenta of the final-state particles and 
the calculation of the local densities are done in 

SUBROUTINE PHASESPACE (RANDOM , ROOTS , P , S , G , I , NG , SWITCH) 
RANDOM : random numbers, 

ROOTS : total scattering energy (after possible ISR), 

P(I, J) : four-momenta of external particles (after possible ISR), 

I = 1, . . . , 7 : numerates the external particles (7 = photon), 

J = 0, . . . , 3 : four- vector components, 
S(I,K) : invariants of pairs of the external particles (S(l,K) = Sik = 2pip k > 0), 
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I, K = 1, . . . , 7 : numerate the external particles (7 = photon), 
G : local density gi of the channel J, 

I : number of the channel, 

NG : number of the phase-space generator, 

SWITCH : decides between event generation, calculation of the local density, and initial- 
ization of the generator, 

SWITCH = : an error occurred (event will be rejected), 
SWITCH = 1 : generation of momenta and invariants, 
SWITCH = 2 : calculation of the local density of the given channel, 
SWITCH = 3 : checks whether the diagram exists (initialization), 
SWITCH = 4 : resulting flag if the diagram exists 

in subtraction. f. The subroutine PHASESPACE is called several times in the subroutine 
SUBTRACTION where the weights of the different subcontributions are calculated. Three 
different phase-space generators are used: 

NG = 1 : four-particle phase-space generator used for tree-level process e + e~ — > Af (in- 
cluding optionally higher-order LL corrections from ISR) or for the IBA to the 
processes e + e~ — > WW — > Af, 

NG = 2 : four-particle phase-space generator for the 4/ part of 0(a) radiative correc- 
tions, 

NG = 3 : five-particle phase-space generator used for the 4/7 part of the 0(a) radiative 
corrections, for the tree-level process e + e~ — >■ 4/7, or for the IBA to the 
processes e + e~ — > 4/7. 

The subroutine PHASESPACE is decomposed into 15 parts which correspond to the differ- 
ent topologies contributing to the processes e + e~ — > Af and e + e~ — > 4/7. There are 5 
topologies for the processes e + e~ — > Af and 10 topologies for e + e~ — > 4/7, as shown in 
Figure 2. 

The initialization of the channels is done by calling 

SUBROUTINE INITPHASESPACE (NUPS , Q , NG , NP ARTICLE , INPUTONSHELL , SRC , 
SIGNAL, NOCUTS) 

NUPS : parameter v used in the propagator mapping of Eq. (9), 

Q(I) : charges of the external fermions (I — 1, . . . , 6), 

NG : number of the phase-space generator, 

NP ARTICLE : number of final-state particles (NP ARTICLE = 4, 5), 

INPUTONSHELL : if INPUTONSHELL = 1 the phase space is generated for on-shell 

W bosons, 

SRC : include special channels for the subtraction functions of the 4/7 part 

of the radiative corrections for SRC = 1 and NP ARTICLE = 5, 



20 



SIGNAL : restricts to subsets of diagrams as given in Table 4, 

NOCUTS : decides whether cuts are included in the event generation (NOCUTS = 0) 

or not (NOCUTS = 1) 

in subtraction. f. The subroutine INITPHASESPACE for the three phase-space generators 
is called in INITSUBTRACTION in subtraction. f, where the options for the phase-space 
generators and for different branches can be found. 

To initialize the phase-space generators, SEARCHGENERATOR (which is just a second entry 
of PHASESPACE) is called in INITPHASESPACE with SWITCH = 3. If the diagram of the 
tested channel exists, SEARCHGENERATOR yields SWITCH = 4, and the channel is used in 
the phase-space generation. The actual selection of the channels is done by inserting the 
external particles EP1, . . . , EP7 (and the emitter EM and spectator SP for special channels 
related to subtraction functions) and the virtual particles IP1, . . . , IP4 in all possible 
ways, and testing whether the vertices of the topological diagram exist or not. This is 
done in the functions V3 and V4 for the three- and four-particle vertices, respectively. 
The vertices of the SM are defined in the functions V3G and V4G. Some additional criteria 
are added to select subsets of diagrams (CHECKSIGNAL), and to avoid double-counting of 
equivalent channels, i.e. channels which have the same propagator mappings. To speed 
up the initialization, the virtual particles are inserted one after another, and as soon as 
a vertex does not exist the diagram is discarded. All informations of the phase-space 
generators are stored in the common block CHANNELS in subtraction. f. 

For SWITCH = 1 the event is generated, i.e. the momenta of the final-state particles are 
calculated from the random numbers RANDOM, the incoming momenta P (1:2, 0:3), and 
the scattering energy ROOTS. For SWITCH = 2 the local density G of the channel I is 
calculated from momenta P, invariants S, and from ROOTS. The invariants S are calculated 
in the subroutine SDET once per event for all local densities. The actual calculations are 
done by calling the subroutines SPROP, DECAY, and PROCESS in subtraction. f: 

SUBROUTINE SPROP (RANDOM ,S,G, MASS , WIDTH , NU , SMIN , SMAX , ONSHELL , SWITCH) 
SUBROUTINE DECAY (RANDOM , Q , PI , P2 , G , SO , SI , S2 , SWITCH) 

SUBROUTINE PROCESS (RANDOM , Ql , 02 , PI , P2 , QT , G , MASS , WIDTH , NU , SO , SI , S2 , T , Tl , 
T2,CCUT, SWITCH) 



RANDOM : random numbers, 

S : invariant generated by SPROP, 

G : local density of the channel, 

MASS, WIDTH : mass and width of the Breit-Wigner propagator mapping of Eq. (7), 

NU : parameter v of the propagator mapping of Eq. (9), 

SMIN, SMAX : minimal and maximal value allowed for the invariant S in SPROP, 

ONSHELL : on-shell approximation for W bosons used in SPROP if ONSHELL = 1, 

CCUT : cut on cosine of scattering angle applied only for 2 — > 2 scattering pro- 

cesses (PROCESS), 
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SUBROUTINE PROCESS: 



SUBROUTINE DECAY: 



SI 



PI = Pi 



01 = qi 



pi 



SO= (P2+P2) 2 
S2 = 



Pi 



T= (qi -Pif 



QT = qi - p 1 



02 = q 2 



Tl = q\ 



T2 = g 2 



S1=P? 



P2 = p 2 PI = p! 




S2=p 2 



P2 =p 2 



SO = g 2 



Fig. 4. The generic diagrams of the 2 — > 2 scattering process and 1 — ► 2 decay and the arguments 
of the corresponding subroutines PROCESS and DECAY in subtraction. f 



SWITCH : see subroutine PHASESPACE, 

the remaining arguments are given in Figure 4. 

The propagator mappings of the massless propagators of Eqs. (7) and (9) and the calcu- 
lation of the corresponding Jacobians are done by calling the functions H and JACOBIAN, 
respectively, in the subroutines SPROP and PROCESS in subtraction. f. For the boosts 
and rotations of momenta the subroutines BOOST and ROTATION in subtraction. f are 
used. 

For the subtraction functions corresponding to the bremsstrahlung process e + e~ — > 4/7, 
special channels are included in the calculation to improve the numerical stability. Since 
these involve the Born cross section to e + e~ — > 4/, a mapping of the 4/7 phase space 
into the 4/ phase space is needed. The mapping from the five- to the four-particle phase 
space and back is performed by calling the subroutines MAP54 and MAP45 in PHASESPACE, 
respectively. The propagators of the Born process e + e~ — > 4f are then smoothened by 
the four-particle phase-space generator. Additional mappings implemented in MAP54 and 
MAP45 account for the smoothening of the remaining peaking behaviour of the subtraction 
functions, i.e. of the splitting functions. 



3.2 The matrix- element evaluation 



In this section, we describe the matrix-element evaluation in the three different modes 
of RacoonWW. Furthermore, we briefly sketch the treatment of higher-order ISR and 
QCD corrections, and the implementation of anomalous gauge-boson couplings. 
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Table 2 

The classification of the four-fermion processes e + e~ — > 4/. Here, / and F denote different 
fermions (/ 7^ F) that are neither electrons nor electron neutrinos (f,F 7^ e~,u e ), and their 
weak-isospin partners are denoted by /' and F' , respectively. The parameters SYM, NC, and CC 
are determined in the subroutine INIT_EPEM in ee4f a_amps . f . 

3.2.1 The tree-level mode (SRC = 0^ 

In the tree-level mode (SRC = 0) RacoonWW calculates either the processes e + e~ — > 4/ 
(SB0RN4 = 1, 2, 3) or the radiative processes e + e" — > 4/7 (SB0RN5 = 1) at lowest order for 
all possible four-fermion final states. The flags SB0RN4 and SB0RN5 are set in the input file. 
The Coulomb singularity to the tree-level process e + e~ — > 4/ can be included by setting 
SCOULTREE = 1, which corresponds to the variant 1 of Section 3.2.2. The treatment of the 
Coulomb singularity in e + e~ — > 4/7 processes is described in Section 3.2.2. 

The different process classes are characterized according to the appearing intermediate 
vector bosons, i.e. whether the reactions proceed via charged- current (CC), or neutral- 
current (NC) interactions, or via both interaction types (see Table 2). The classification 
can be performed by considering the quantum numbers of the final-state fermion pairs. 
As shown in Table 2, we identified 11 classes of processes. This classification does not 
change when a photon is radiated. In the following, we often use the shorthands CC11 
and CC03. CC11 denotes the set of up to 11 Feynman diagrams contributing to the 
simplest CC process in Table 2 (SYM = 0), which is the smallest gauge-invariant subset 
of diagrams including W-pair production. CC03 denotes the 3 diagrams that include 2 
resonant W bosons, which is a non-gauge-invariant subset of CC11 (see also Ref. [10]). 

The calculation of the matrix elements to the tree-level processes e + e" — > 4/ and e + e~ — > 
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b) Non-abelian graph 
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fc(Pc, <?c) 
fdiPd, cr d ) 

fe(Pe,(Te) 
ffiPf^f) 



Fig. 5. Generic diagrams for e + e — ► 4/ 



4/7 has been described in detail in Ref. [4]. For e + e~ — > 4/ the matrix elements of all 
Feynman diagrams can been reduced to two generic functions that are related to the two 
basic graphs shown in Figure 5. For e + e~ — > 4/, the calculation is similar to the one in 
Ref. [27]. For e + e~ — > 4/7, all graphs that are obtained from one of the generic graphs 
in Figure 5 by attaching one additional real photon in all possible ways are combined 
into one generic function. Finally, also all amplitudes to e + e~ — > 4/7 can be constructed 
from only two generic functions. For the helicity amplitudes corresponding to the generic 
contributions, concise results have been obtained by using the Weyl-van der Waerden 
(WvdW) formalism (see Ref. [28] and references therein). 

The helicity amplitudes to e + e~ — > 4/ are calculated in the function M2_EPEM_4F by calling 
either the subroutine M_CC or M_NC or both, depending on the process under consideration: 

SUBROUTINE M_CC (M , SGN , IF1 , IF2 , IF3 , IF4 , IF5 , IF6) 
SUBROUTINE M_NC (M , SGN , IF1 , IF2 , IF3 , IF4 , IF5 , IF6) 
M(-l:l,. . . ,-1:1) : helicity amplitudes, 
SGN : relative sign of different diagrams, 



IF1,. ..,IF6 



: specify the final-state fermions. 



The helicity amplitudes are expressed in terms of WvdW spinor products, which are 
calculated in the subroutine SETPRODS. Note that SETPRODS must to be called before 
the matrix elements are evaluated by the subroutines M2_EPEM_4F or M2_EPEM_4FA. The 
helicity amplitudes for a four-fermion final state, M, are built up by using two generic 
functions related to the graphs of Figure 5 as calculated in the subroutines ADDMVV for 
the abelian diagram and ADDMVVV for the non-abelian diagram of Figure 5, respectively, 
with the appropriate choices for the parameters SGN, IF1,. . . ,IF6, and by specifying the 
intermediate gauge bosons V, Vi, V2 of Figure 5 by the variables IV, I VI, IV2, which are 
arguments of ADDMVV and ADDMVVV. 

The generic functions for e + e~ — > 4/7 can be constructed in a similar way. The helicity 
amplitudes to e + e~ — > 4/7 are built up in the subroutines M_CCA and M_NCA. Again 
the contributions of all diagrams can be taken into account by calling the corresponding 
subroutines ADDMAVV and ADDMAVVV and making the appropriate choices for the parameters 
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SGN, IF1,...,IF6, IV, IV1, IV2. 



For the NC (SYM = 0, 1) and CC/NC (SYM = 0) type of processes with pure hadronic final 
states, one of the gauge bosons in the diagram of Figure 5 can be a gluon. In RacoonWW 
these contributions are taken into account when choosing SQCDEPEM = 1 in the input file 
for both processes e + e~ — > 4/ and e + e~ — > 4/7. Then, the complete matrix elements for 
the fully hadronic channels result from the sum of the purely electroweak and the gluon- 
exchange contributions. The contribution of the gluon-exchange diagrams are calculated 
in the subroutines MG_NC and MG_NCA for the process classes e + e~ — > 4/ and e + e~ — > 4/7, 
respectively. 



3.2.2 The improved- Born approximation (IBA) mode (SRC = 2,3,) 

RacoonWW provides an IBA for the process e + e" — > WW — > 4/ (SRC = 2) including 
only diagrams of the CC03 class, and to the radiative processes e + e~ — ■> 4/7 (SRC = 3), as 
described in Refs. [7] and [8], respectively. For both process classes universal leading elec- 
troweak corrections are taken into account. Universal radiative corrections are corrections 
that are either connected to specific important subprocesses, such as collinear photon 
emission, or that originate from large corrections related to renormalization effects, like 
running couplings, and involve characteristic enhancement factors. Owing to their univer- 
sality such corrections can be related in a simple way to the lowest-order matrix element 
of the underlying process. 

Near the W-pair-production threshold, a particularly important correction is induced by 
the Coulomb singularity. This arises from diagrams where a soft photon is exchanged 
between two nearly on-shell W bosons close to their kinematical production threshold 
and results in a simple factor that depends on the momenta k± of the W ± bosons [17,29], 



do- C oui= 1 + 5coui(s, k.,k_)g(0) 



da 



Born? 



r 1 ,2 ,2\ t 

dcouifs, k + , k_) = — =- Im 



P 



In 



'P-P + A 



M 



P + P + A 



M 



y/s 2 + k\ + ki- 2sk\ - 2sk 2 _ - 2k\kl 



^^^-iMw^ Am 



\kl - kl 



(24) 



with the fine-structure constant a(0), and s = (k + + k^) 2 . The auxiliary function 



-„\ 2 



l(P) = (l~p 2 ) 



(25) 



is introduced to restrict the impact of <5c ou i to the threshold region. 

A further ingredient in the IBA is the LL contribution induced by ISR. This is included by 
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using the structure-function approach [18], where the full IBA matrix elements squared 
to e + e~ — > WW — > 4/ and e + e~ — > 4/7 are multiplied by the ISR structure functions. 
Their scale dependence can be used to adjust the IBA to the full correction, but also to 
estimate the intrinsic uncertainty of the IBA by choosing different values for Q. 

When the IBA to e + e~ — > WW — > 4/ is calculated, the (off-shell) improved Born matrix 
element squared described by 



FUNCTION M2_B0RN_0FFSH (P , QCC03 , QGZ) 
P(l:7,0:3) : external momenta (after possible ISR), 

QCC03 : chooses between the lowest-order CC03 cross section (QCC03=0) and the 

IBA for the CC03 cross section (QCC03=1), 

QGZ : decides whether the Z-boson width, GZ, is used in the Z propagator 

(QGZ = 1) or not (QGZ = 0) 

in eeWW4f _amps . f is convoluted with the ISR structure functions. The Born matrix ele- 
ment includes all contributions from the CC03 class of diagrams (QCC03 = 1) which have 
been modified in such a way that the universal renormalization effects induced by the 
running of a and by Ap are absorbed. This is achieved [30] by the replacements 

e 2 

— -> AV2G^, e 2 -> Ana(s) (26) 



in the amplitudes, which implies that weak-isospin exchange involves the coupling G^M^ 
and pure photon exchange the coupling a(s). The running of the electromagnetic coupling 
is induced by light (massless) charged fermions only, i.e. we evaluate a(s) by 

aW = gCjg) (27) 



where a(M|) = ALPHAZ is set in the subroutine PARAMETER in public . f . Furthermore, the 
Coulomb singularity is taken into account, and naive QCD corrections can be optionally 
included. The W-boson width is treated in the constant-width scheme. The complex Z- 
boson mass is used (QGZ = 1) in the Z propagators in order to regularize the Z resonance 
below the W-pair production threshold, which may be reached after ISR. 

The IBA for e + e~ — > 4/7 processes is based on the full matrix elements of the tree-level 
mode in the G M scheme. Moreover, leading ISR corrections are included via structure 
functions, and the naive QCD correction can be optionally included. Again, the constant- 
width scheme is used. As the tree-level mode, the IBA is not restricted to a certain class 
of processes, and the subset of diagrams included is governed by the choice of the flags 
SSIGEPEM5 and SQCDEPEM in the input file. For e + e" — > 4/7 the naive implementation 
of the Coulomb singularity involves some freedom, since the momenta of the W bosons 
depend on the photon momentum in different ways for photon radiation off initial-state 
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and final-state particles. To control this ambiguity the Coulomb factor can be evaluated 
in two different ways: 

(f) In the first variant (SCOULTREE = f), the complete matrix element squared is multi- 
plied with the Coulomb correction factor using k + = k 3 + k 4 and k_ = k 5 + k 6 . In this 
way the correct Coulomb correction is multiplied to all diagrams where the photon 
is radiated off initial-state particles. 

(2) In the second variant (SCOULTREE = 2) this prescription is improved by choosing 
the correct Coulomb factor to all contributions involving the most resonant W + and 
W~ bosons as determined from invariant masses in the final state. To this end, the 
W-boson momenta entering the Coulomb correction factor are fixed as 



(k 3 + fc 4 , k 5 + fc 6 ) 


for 


A 34 


< A347, A 56 


< A 567 , 






(k s + k 4 + k 7 , k 5 + k 6 ) 


for 


A 34 


> A347, A 56 


< A 567 
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> A347, A 56 


> A 56 7, 


A34 7 


< A 567 , 


(k 3 + k 4 , k 5 + k 6 + k 7 ) 


for 


A 34 


< A 347 , A 56 


> A 567 


or 








A 34 


> A 347 , A 56 


> A 567 , 


A3 47 


> A 567 , 



(28) 

where = \{ki + kj) 2 — M^\ and A^-j = \{ki + kj + ki) 2 — M^\. Thus, the correct 
Coulomb correction factor is applied to all dominating doubly-resonant contributions. 

Finally, for the IBA the W-boson width rv is calculated in lowest order using the 
G M scheme. This choice guarantees that the "effective branching ratios", which result 
after integrating out the decay parts, add up to one when summing over all channels. Of 
course, if naive QCD corrections are taken into account by multiplying with (l + a s /ir) for 
each hadronically decaying W boson, for consistency these QCD factors are also included 
in the calculation of the total W-boson width. 

3.2.3 The radiative correction (RC) mode (SRC = 1) 

In the RC mode (SRC = I) the complete electroweak corrections to e + e~ — > WW — > 4/ 
in DPA are calculated as described in detail in Ref. [6]. Including 0(a) corrections, the 
cross section to e + e~ — > 4/(7) is composed of the following contributions 

d<r = datfr 41 + d^r^ + d<7 e+e "^. (29) 

Here da^ n ^ and do- c+c ~ >4 ^ 7 are the lowest-order cross sections to e + e~ — > 4/ and to 
e + e~ — > 4/7, respectively, and da^ t denotes the virtual one-loop corrections. 

The e + e~ — > 4/7 process (SB0RN5 = 1) represents the bremsstrahlung contribution, where 
the soft and collinear singularities are either subtracted from the matrix element squared 
(subtraction branch) or avoided by imposing technical cuts on the photon phase space 
(slicing branch). In the RC mode the calculation of the e + e~ — > 4/7 process is always 



27 



restricted to the CC11 class of diagrams (SSIGEPEM5 = 5). Simultaneously, the tree-level 
process e + e~ — > 4/ (SB0RN4 = 1) can be calculated including all diagrams or only a subset 
of diagrams depending on the choice for SSIGEPEM4 in the input file. 

Both the virtual and the real corrections involve soft and collinear singularities. These 
singularities are extracted by separating the cross sections into finite and singular parts: 



i e+e~^4/ i e+c — >4/ . i c+c — >4/ 

UC'virt ^^virt, finite ' ^^virt,sing i 

da c+c ~^ = dalX^ 1 + d^"^ 7 - (30) 
In RacoonWW, the DPA is only applied to the finite part of the virtual corrections, 



i e+e-^4/ i e+e-^WW^4/ _ , e+e^WW-»4/ _ i c+c^WW^4/ /oi\ 

a<T virt, finite — ¥ a °virt,finite,DPA — a0 virt,DPA a0 virt,sing,DPA ■ 

The doubly-resonant virtual corrections are composed as 



d^ D PA WW - 4/ = ^d$ 1/ ^Ro 



Born, DPA j UJ v 1 virt,fact,DPA 



+ 



^'Born^PA 



"nfact.DPA 

(32) 



where fa^DPA^ 4 ^ denotes the matrix element for the factorizable virtual corrections 

and ^nfact dpa i s t ne factor describing the non-factorizable virtual corrections. The matrix 
element squared to both contributions are calculated in eeWW4f _DPA . f by 

FUNCTION M2(P,P0N) 

P ( 1 : 7 , : 3) : four-momenta of external particles, 

P0N(1:7,0:3) : four-momenta of external particles for on-shell W bosons. 

The full four-fermion phase space (with off-shell W bosons) is generated, but the matrix 
elements for the DPA part are evaluated with kinematics for on-shell W bosons. This 
necessitates a mapping of the momenta of the final-state fermions such that the momenta 
of the W bosons become on shell. This is done in eeWW4f_DPA.f by 

SUBROUTINE 0FT00N(P,P0N,MW2,N) 

P,P0N : see function M2, 

MW2 : W-boson mass squared, 

N = 1, . . . ,4 : different versions of the on-shell mappings (see Section 4.1). 

The factorizable contribution to the matrix element squared is calculated in the function 
DM2 (F) in eeWW4f _amps . f in terms of invariant coefficients F(l:10,-l:l) which multiply 
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so-called standard matrix elements, M.°, asQ 



n=l 

(33) 

where k± are the on-shell momenta of the W ± bosons, and t = (p + — k + ) 2 . It is important 
to realize that the invariant coefficients F°(s,t) = F(n,a) depend only on the scalar 
products of the momenta of the on-shell-projected W bosons, but not on the momenta 
of their decay products. For the default choice QREAL = 1, the real part of F is taken, 
i.e. only the real part of the factorizable virtual corrections is included. The form factors 
to e + e" -> WW -> 4/ are calculated in the subroutine CALCFORM for QFAST = and in 
FASTFORM for QFAST = 1 in eeWW4f _DPA . f . The standard matrix elements are calculated 
in the subroutine CALCSME in eeWW4f _amps . f . 

The actual numerical evaluation of the invariant coefficients F£(s,t) can be rather CPU- 
time consuming, since the expressions for the one-loop contributions to these coefficients 
are rather involved (see Refs. [31,32]). Moreover, the employed Passarino-Veltman reduc- 
tion [33] of tensor integrals to scalar integrals breaks down at the boundary of the phase 
space, necessarily leading to numerical instabilities in the very forward and backward 
directions of the W-production angle 9. In order to solve these problems of CPU time 
and numerical instabilities, Racoon WW uses as default (QFAST = 1) a numerical re- 
construction of the coefficients by making use of an expansion in Legendre polynomials 



E 

1=0 



id" 
t n 



[s)Pi (cos I 



(34) 



which involves only simple algebra. Thus, its evaluation is extremely fast. In practice, 
relatively few generalized Fourier coefficients c^ ; (s) are needed; for example, taking I = 
0, . . . , 20 reproduces the full calculation of F°(s,t) within roughly six digits for moder- 
ate scattering angles and LEP2 energies. The generalized Fourier coefficients are calcu- 
lated via Gaussian integration during the initialization in the subroutine CFOURIER in 
eeWW4f _DPA. f . Note also that the generalized Fourier series remains stable in the forward 
and backward directions where the original evaluation of F£(s,t) breaks down. The ex- 
plicit factor t in Eq. (34) was included in order to account for a t-channel pole in some 
of the F°(s,t); without this factor the expansion is less efficient and requires much more 
terms in the sum over I in Eq. (34). The Legendre polynomials Pi(cos6) are calculated in 
eeWW4f_DPA.f by the subroutine LEGP. 



1 In Ref. [6] the sum over n runs from 1 to 7, since the terms for n = 8, 9, 10 could be expressed 
in terms of the first 7 terms. However, in RacoonWW the version with 10 invariant coefficients 
is used. 
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As described in Ref. [6] (and references therein), the virtual non-factorizable corrections 
can be written in terms of a correction factor 5^act dpa to the lowest-order cross section. 
This factor is calculated in eeWW4f_DPA.f by the function DNONF with the possibility of 
choosing between different implementations (QNF = 1,2,3), as described in Section 4.1. 

The remaining part of the electroweak radiative corrections, i.e. the soft and collinear 
parts of virtual and real photon radiation and the hard bremsstrahlung contribution are 
calculated exactly, i.e. using off-shell kinematics and based on the full matrix elements to 
e+e~ -> 4/ and e + e~ -> 4/7 (SSIGEPEM5 = 5) of the tree-level mode. The matching of the 
virtual photon corrections calculated in DPA with the exact real photon radiation cannot 
be done unambiguously as far as the IR-finite terms are concerned. The definition of the 
virtual singular cross section is a matter of convention, since finite terms can be redis- 
tributed between singular and finite parts. This freedom, which is beyond DPA accuracy, 
is reflected by the different choices provided by Racoon WW in defining da^ fi ^ t ^p A 4 ^ 
(QSOFT = 1, 2, 3 and, for the subtraction branch, SSUB = 1, 2). 

The singular part of the virtual corrections, dcr^ si ^ g 4 ^, factorizes into the lowest-order 
cross section and a simple correction factor. The correction factor of the part which is 
subtracted, si^DPA* 4 ^ is calculated in the slicing branch by the function DS0FT1 in 
eeWW4f_DPA.f for QSOFT = 1,2 or by calling SOFTFAC in kern.f for QSOFT = 3. In the 
subtraction branch, the part which is subtracted is calculated by the subroutine SUB4 
in subtraction. f for QSOFT = 1,2,3. The corresponding lowest-order cross section to 
e + e" — > WW — > 4/ is calculated in the subroutine CALCSME in eeWW4f _amps . f . 

The soft and collinear singularities in the real part, d<7g ing ^ 4 ^ 7 5 can be split off, and 
the integration over the singular regions of the photon phase space can be performed 
analytically. The result can be written as a convolution of a structure function with the 
lowest-order cross section dcrp^ - * . When adding dal ir l s ^^ and d<Tg in g ~^ 4 ^ 7 , all soft 
singularities and all collinear singularities associated with the final state cancel, and the 
resulting "singular" cross section, 

1 e+e~^4/ 1 c+c~— >4/ . j c+c~-^4/7 ("IK} 

uo virt+real,sing uo virt,sing ' uo sing J yoo J 



contains, apart from finite terms, only collinear singularities associated with the initial 
state, i.e. leading logarithms of the form ln(s/rrig), at least for sufficiently inclusive ob- 
servables. 

Finally, we arrive at the master formula for the cross section: 

OCT — QCr Bom -|- U(7 virt)finite)D p A -+- U(7 virt+real)Sing "I" Uff finite . {OU) 



Since the contribution dcr^ + ^^ si is not treated in DPA, the LL photonic corrections 
resulting from ISR are treated exactly in our approach. 
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As a result, in the RC mode Racoon WW calculates several weights corresponding to 
processes with 4/ and 4/7 kinematics. The cross sections da^ n ~* 4 ^ , dal ir l ^^dpa^ 1 

dCyirt+reaifsing' and d<7g^ e of Eq. (36) are obtained from the corresponding 4/ weights 
and the 4/7 weight, respectively. 



3.2.4 Higher-order leading-logarithmic (LL) corrections 

Higher-order LL corrections from ISR can be taken into account in all three modes by 
setting SISR = 1 in the input file and are implemented in LL approximation via the 
structure-function approach (see Section 5 of Ref. [6]), i.e. by convoluting the tree- level or 
IBA cross section with the structure functions of Refs. [9,18]. The ISR structure functions 
depend on the QED splitting scale, Q = SCALE, which is not fixed in LL approximation and 
has to be set to a typical momentum scale of the process. In the RC mode 0(a) corrections 
from ISR are calculated completely, and the 0(a) LL corrections are implicitly included. 
In order to avoid double counting of radiative corrections, the LL 0(a) contributions from 
structure functions multiplied by the CC11 Born cross section are subtracted from the 
full 0(a) corrections. Thus, for SRC = 1 and SISR = 1, the dependence on the QED scale, 
Q, only enters in 0(a 2 ), i.e. at the two-loop level. 

In the slicing branch of RacoonWW the ISR structure functions are calculated in kern . f 
by 

FUNCTION STRUCI (XX , DXX , LLE , BET , ALP) 
XX,DXX : random numbers, DXX = 1 - XX, 

LLE,BET,ALP : BET = ALP (LLE - 1)/tt, LLE = 2 ln(SCALE/m e ), ALP = a(0). 

In the subtraction branch they are obtained by calling 

SUBROUTINE ISR (RANDOM , XI , X2 , STFCT , WISR , SWITCH) 
RAND0M(1:2) : random numbers, 

XI, X2 : energy fraction of incoming particles (after ISR), 

STFCT : value of the structure functions only, 

WISR : total weight for generation of XI, X2 including structure functions, 

SWITCH : auxiliary flag (SWITCH = if an error occurs), 

in subtraction.!. By default, SCALE is set to the CM energy, SCALE = y/s, in the 
subroutine INITIALIZE in racoonww.f. 



3.2.5 Naive QCD corrections 

In all three modes of RacoonWW, the so-called naive QCD corrections can be included 
by setting QQCD = 2 in the input file. QCD corrections enter the processes with hadronic 
final-states at two different places. On the one hand, the hadronic W-boson width receives 
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a QCD correction, namely a factor (1 + a s /7r) [34,35], if the W-boson width is calculated 
(QGW = 4). This affects the total W-boson width T-w in the W propagators. On the other 
hand, each hadronically decaying W boson receives a QCD correction to the Wqq' vertex. 
If the full phase space for gluon emission is integrated over, this correction reduces to 
a multiplicative correction factor (1 + a s /ir) for each hadronically decaying W boson in 
0(a s ). The application of this inclusive factor to distributions is usually called naive 
QCD correction. In the total cross section for W-pair production naive QCD-correction 
factors cancel against the corresponding factors in the W-boson width in the resonant 
W propagators. 



3.2.6 Anomalous triple and quartic gauge-boson couplings 

In RacoonWW anomalous triple gauge-boson couplings (ATGCs) to e + e~ — > 4/ and 
anomalous quartic gauge-boson couplings (AQGCs) to e + e~ — > 4/7 have been imple- 
mented. More precisely, the non-standard couplings are introduced in the lowest-order 
matrix elements, which can be dressed by initial-state radiation, while the contributions 
from the remaining radiative corrections add linearly to the non-standard effects. The 
ATGC contributions are taken into account when SATGC = 1 is set in the input file 
(see also Ref. [15]). The charged ATGCs are defined as in Ref. [36]. The conventions for 
the neutral ATGCs follow those of Refs. [37,38]. The ATGC contributions to the ma- 
trix elements are expressed in terms of five generic functions corresponding to the ZVV 
and VWW couplings with V — 7, Z, which have been implemented in the subroutines 
ADDMATGCVVZ and ADDMATGCVWW of ee4f a_amps . f , respectively. These ATGC contribu- 
tions are added to the generic CC and NC helicity amplitudes in the subroutines M_NC 
and M_CC. When SATGC = 1 is chosen the result depends on the ATGC parameters Ag\ , 
Ak v , \ v , k v , \ v , g± 5 , fY 5 , h\ 3 with V = 7, Z, which are specified in the input file (see 
Section 4.1). 

As described in detail in Ref. [8], AQGCs that involve at least one photon have been 
implemented in RacoonWW. They are taken into account when QAQGC = 1 is set in 
the input file. Only genuine AQGCs are considered, i.e. all operators that contribute also 
to ATGCs, such as the quadrilinear part of the well-known operator F^W+' P W~^, are 
omitted. Imposing in addition custodial SU(2) C invariance [39] to keep the p parameter 
close to 1, only operators of dimension 6 or higher occur. Following Refs. [8,40,41,42,43] 
dimension-6 operators for genuine AQGCs are considered that respect local U(l) em in- 
variance and global custodial SU(2) C invariance. These symmetries reduce the set of such 
operators to a phenomenologically accessible basis. More general AQGCs were discussed 
in Ref. [44]. The AQGC contributions are expressed in terms of three generic functions 
which correspond to the r j r yVV and 7ZWW couplings with V = W, Z, respectively. They 
are added to the generic Standard Model functions in M_CC and M_NC by calling the sub- 
routines ADDMAAQGCAAWW, ADDMAAQGCAAZZ, and ADDMAAQGCAZWW in ee4f a_amps . f . When 
QAQGC = 1, the result depends on the AQGC parameters a 0iC>n /A 2 and a , n /A 2 specified 
in the input file (see Section 4.1). 
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1 
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* 


* 
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racoonww.f 


on-shell projection 
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* 


* 


* 
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QSOFT 


racoonww.f 


photon DPA part 


* 


* 


* 


* 


1-3 


OPT 


slicing.f 


optimization 


0,1 


0,1 


0,1 


0,1 


0,1 


SSUB 


subtraction. f 


subtr. functions 


* 


* 


* 


* 


1,2 



Table 3 

The allowed choices for the options in the different modes of RacoonWW 
4 The usage of RacoonWW 



In the sequel we describe the input parameters and flags for the options of RacoonWW, 
starting with those that are set in the input file.[^]A list of all possible values of the flags of 
RacoonWW discussed below is given in Table 3. Some more information on the options 
and details about the format in the input file can be found in the file README. 

2 The flags SB0RNG5, SSIGEPEMG5, SREC0MBG, and the parameters PRECDMBG are not listed here 
and are not discussed any further, since they only affect the QCD corrections beyond the naive 
approximation, which are not supported by RacoonWW yet. 
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4-1 Input and initialization of RACOON WW via the input file 

OUTPUTFILE: a character string that specifies the name of the output file. 

ENERGY: the total CM energy in GeV (ENERGY = y/s). 

NEVENTSW : number of weighted events. 

The number of weighted events must be at least NEVENTSW = 2 x 10 6 in the tree-level and 
IBA modes and NEVENTSW = 10 7 in the RC mode to guarantee that the multi-channel 
integration yields reliable estimates for the MC integration error. When unweighted 
events are generated (NEVENTSUNW > 0), NEVENTSW is the number of weighted events 
used to initialize the unweighting procedure; these weighted events are used to determine 
the maximal weights (see Section 4.4). 

SMC : flag for the slicing and subtraction branches of RACOONWW: 
SMC = 1 : chooses the slicing branch in its default setup, 
SMC = 2 : chooses the subtraction branch, 

SMC = 3 : chooses a slightly modified version of the slicing branch, which has been 

introduced to minimize the number of negative weights in the RC mode. 
In the RC mode, RACOONWW provides an unweighting procedure only for SMC = 3. 
More details can be found in Section 4.4. 

SB0RN4 : flag for the tree-level e + e~ — > 4/ cross section: 

SB0RN4 = : the tree-level contribution to e + e~ — > 4/ is not included, 
SB0RN4 = 1 : calculates the off-shell Born e + e^ — > 4/ cross section, 

SB0RN4 = 2 : calculates the cross section to e + e~ — > WW — > 4/ in DPA with phase- 
space for off-shell W bosons, 

SB0RN4 = 3 : same as SB0RN4 = 2 but with phase-space for on-shell W bosons. 
For SB0RN4 = 1 the selection of a subset of diagrams is governed by the flag SSIGEPEM4 
(see Table 4). For SB0RN4 = 2, 3 the Born cross section is based on the CC03 diagrams 
only. While for SB0RN4 = 2 the full phase space for off-shell W bosons is generated, the 
matrix element is calculated with on-shell momenta which result from a projection of 
the phase space with off-shell W bosons onto a phase space with on-shell W bosons. 
For SB0RN4 = 3 both the phase space and the matrix element squared are calculated 
with on-shell kinematics. In both cases the Z-boson width is set to zero. The flags 
SB0RN4 = 2,3 should only be used for checks and comparisons but not for the generation 
of physical results. Moreover, SB0RN4 = 2, 3 is only a valid choice for CC processes or 
mixed processes with SYM = 0, 3 (see Table 2). 

SB0RN5 : flag for the tree-level cross section to e + e~ — > 4/7: 

SB0RN5 = : the tree-level contribution to e + e~ — > 4/7 is not included, 

SB0RN5 = 1 : calculates the tree-level cross section to e + e~ — > 4/7. 
In the tree-level mode with SB0RN5 = 1, Racoon WW calculates the Born contribution 
to e + e" — > 4/7. In the RC mode (SRC = 1) and in the IBA mode to e + e" — > 4/7 
(SRC = 3), RACOON WW always sets SB0RN5 = 1 in racoonww.f independent of the 
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choice in the input file. In the slicing branch (SMC = 1,3), when radiative corrections 
to e + e~ — > WW — > 4/ are calculated (SRC = 1), SB0RN5 = 1 calculates the hard- 
photon bremsstrahlung contribution, i.e. photon radiation away from soft and collinear 
singularities. 

SISR: flag for the inclusion of higher-order ISR: 

SISR = : higher-order LL corrections from ISR are not included, 

SISR = 1 : higher-order LL corrections from ISR are included in the LL approxima- 
tion up to order 0(a 3 ) for SMC = 1,2 and up to order 0(a 2 ) for SMC = 3. 
Higher-order LL corrections from ISR are taken into account by convoluting the tree- 
level cross section with the structure functions of Refs. [9,18] (see also Section 3.2). 
In the RC mode the ISR beyond 0(a) is switched on/off by the choice of SISR, i.e. 
0(a) ISR contributions are always included as part of the 0(a) corrections even when 
SISR = 0. In the RC mode with SISR = the full 0(a) ISR contribution is included for 
all diagrams of the CC11 process class. For SRC = 1 and SISR = 1 the LL contribution 
to the ISR is taken into account for the full four-fermion process, i.e. also for the back- 
ground diagrams, depending on the choices for SSIGEPEM4. When ATGCs are considered 
(SATGC =1), the complete ISR corrections in LL approximation (0(a)+higher-orders) 
to the additional contributions involving anomalous couplings is also switched on/off by 
the choice of SISR. When calculating the IBA (SRC = 2,3), or the tree-level processes 
e+e" -> 4/ or e + e" -> 4/7 (SRC = 0, SISR = 1, SB0RN4 = 1 or SB0RN5 = 1), the full 
corrections from ISR are taken into account in LL approximation up to order 0(a 3 ). 

SRC : flag for choosing the tree-level, IBA, or RC mode: 

SRC = (tree-level mode) : calculation of the tree-level processes e + e~ — > 4/ 

(SB0RN4 = 1) or e+e" -> 4/7 (SB0RN5 = 1), 

SRC = 1 (RC mode) : 0(a) radiative corrections to e + e~ — > WW — > 4/ are 

included in DPA, 

SRC = 2 (IBA mode) : IBA for e+e" -> WW -> 4/ is calculated for the CC03 

class of diagrams, 

SRC = 3 (IBA mode) : IBA for e + e~ — > 4/7 is calculated based on the full ma- 
trix element. 

SCOULTREE: flag for the inclusion of the Coulomb singularity in the tree-level and IBA 
modes: 

SCOULTREE = : the Coulomb singularity to e + e~ -> 4/ and e + e" -> 4/7 is not 
included, 

SCOULTREE = 1 : the off-shell Coulomb singularity to e + e~ — > 4/ in the tree-level 
mode and to e + e~ — > 4/7 in the tree-level and IBA modes is in- 
cluded, 

SCOULTREE = 2 : the off-shell Coulomb singularity to e + e" -> 4/7 in the tree-level 

and IBA modes is included. 
For SCOULTREE = 1 the Coulomb singularity is calculated with invariant masses derived 
from the four-momenta of the final-state fermions only. For SCOULTREE = 2, which is 
only a valid choice for the processes e + e~ — > 4/7, the Coulomb singularity is calcu- 
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lated with the invariant masses of those W + and W~ bosons which are most resonant. 
Note that the invariant masses of the W bosons may or may not include the photon 
momentum depending whether the photon is emitted from the initial or final state (see 
Section 3.2). The choices SCOULTREE > are only valid for WW mediated processes, 
i.e. for CC processes and mixed processes with SYM = 0, 3. 

QNF : flag for the inclusion for the Coulomb singularity in the RC mode: 
QNF = 1 : DPA Coulomb singularity, 

QNF = 2 : off-shell Coulomb singularity with DPA Born cross section, 

QNF = 3 : off-shell Coulomb singularity with (off-shell) CC03 Born cross section. 
This flag is only effective in the RC mode. According to the choice of QNF, different im- 
plementations of the Coulomb singularity are included in the virtual non-factorizable 
corrections (the real non-factorizable corrections are implicitly contained in the full 
matrix-element calculation of e + e~ — > 4/7). The off-shell Coulomb singularity is im- 
plemented as defined in Ref. [17]. In Ref. [6] different versions of the DPA have been 
investigated by modifying the implementation of the DPA. Different treatments of the 
Coulomb singularity, QNF = 2 ("def") and QNF = 1 ("Com"), are one of the considered 
options. The recommended setting is QNF = 3. 

QREAL : flag for the inclusion of the imaginary parts of the virtual electroweak corrections: 
QREAL = : imaginary parts of the one-loop corrections are included, 

QREAL = 1 : imaginary parts of the one-loop corrections are not included. 
This flag is only effective in the RC mode. The imaginary parts of the one-loop correc- 
tions contribute only significantly to distributions that depend on the azimuthal decay 
angle of one of the W bosons. If such distributions are considered, QREAL has to be set 
to QREAL = (imaginary parts are not neglected). 

QALP : flag to specify the input-parameter scheme: 
QALP = : a(0) scheme, 
QALP = 1 : a(M z ) scheme, 
QALP = 2 : Gf, scheme. 

Independent of the choice of the input-parameter scheme, the electroweak corrections 
are always calculated by using the fine structure constant a(0), i.e. the relative (elec- 
troweak) corrections (dcr — da"B rn)/da"Born are always proportional to a(0). Not all 
input parameters defined in the subroutine PARAMETER in public.f are used in the dif- 
ferent input-parameter schemes, e.g. ALPHAZ is only relevant for QALP = 1, GF only for 
QALP = 2. While the small fermion masses have to be non-zero for all input-parameter 
schemes, their actual values matter only if QALP = 0, except for the electron mass that 
is used in the structure functions. 

QGW : flag for the W-boson width: 
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QGW = : the total W-boson width is an input parameter, set in the subroutine 
PARAMETER in public.f, 

QGW = 1 : the tree-level W-boson width is calculated, 

QGW = 2 : the IBA W-boson width is calculated, 

QGW = 3 : the one-loop W-boson width is calculated without 0(a s ) corrections, 
QGW = 4 : the one- loop W-boson width is calculated with 0(a s ) corrections. 
When the one-loop W-boson width is calculated (QGW = 3,4 and SRC = 0,1), the 
full electroweak 0(a) corrections in the chosen input-parameter scheme (controlled 
by the flag QALP) are included. When QGW = 2, the W-boson width is calculated in 
IBA, i.e. the tree-level W-boson width in the G M scheme (QALP = 2) including 0(a s ) 
corrections. The only allowed options for SRC = 2,3 are QGW = 0,2. In racoonww.f 
the partial decay widths for the decays into leptons and quarks are also calculated 
using the same options as for the calculation of the amplitude. In the output file, 
Racoon WW provides the "effective" branching ratios: leptonic width/total width, 
hadronic width/total width, and (leptonic + hadronic width) /total width. The widths 
in the numerators are calculated with the same options as the DPA matrix element 
for e + e~ — * WW — ■> 4f. Since the total W-boson width and the partial widths are not 
necessarily calculated with the same options, the total branching ratio can be different 
from one. The effective branching ratios are not used in RacoonWW and are only 
provided for information. 

QPROP : flag for choosing a width scheme: 

QPROP = : gauge-boson widths are set to zero, 

QPROP = 1 : constant width is used in all propagators, 

QPROP = 2 : constant width is used only in time-like propagators, 

QPROP = 3 : running width is used (no width in space- like propagators), 

QPROP = 4 : constant width is used in all propagators and the weak mixing angle is 

calculated from the complex W- and Z-boson masses. 
For SRC 7^ 0, QPROP = 1 is automatically chosen in the subroutine INITIALIZE in 
racoonww.f. 

SSIGEPEM4: chooses a subset of diagrams for the tree-level processes e + e~ — > 4/ (see 
Table 4). 

SSIGEPEM5 : chooses a subset of diagrams for the tree-level and IBA processes e + e~ — > 4/7 
(see Table 4). 

QQCD : flag for the inclusion of QCD corrections: 
QQCD = : QCD corrections are not included, 
QQCD = 2 : naive QCD corrections are included. 

In RacoonWW only so-called naive QCD corrections can be included yet, i.e. for each 
hadronically decaying W boson Racoon WW applies a global factor (1 + a s /n) to the 
differential cross sections. 

FNAME(3:6): particle names for the specification of the final-state fermions: 
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SSIGEPEM4 
SSIGEPEM5 


subset of diagrams included 


n 


all pi pplrmAfpP! V n l pi crva in 


1 


only electroweak WW signal diagrams (CC03) 


2 


only electroweak ZZ signal diagrams (NC02) 


3 


only electroweak ZZ, Z7, and 77 signal diagrams 


4 


only electroweak WW, ZZ, Z7, and 77 signal diagrams 


5 


only diagrams corresponding to the CC11 process class 



Table 4 

Possible choices for subsets of Feynman diagrams as governed by SSIGEPEM4 and SSIGEPEM5 
FNAME = u, d, s, c, t, b, e, mu, tau, nu_e, nu_mu, nu_tau. 

When specifying the final-state particles the order of Table 2 must be respected. 

PP , PM: degrees of positron- and electron-beam polarization; for unpolarized beams choose 
PP = PM = 0. 

SRECOMB: flag for the photon-recombination scheme: 
SRECOMB = : no recombination is performed, 

SRECOMB = 1 : TH recombination scheme is applied (see also Ref. [6]), 
SRECOMB = 2 : EXP recombination scheme is applied (see also Ref. [6]), 
SRECOMB = 3 : TH recombination scheme is applied including photon recombination 
also with the beams. 

When SRECOMB = 1,2,3 the parameters PRECOMB (1 : 4) must be specified in the input file 
(see below). For SRECOMB = 3 no separation cuts are required and PRECOMB(l) = can 
be chosen. Note that if the unweighting procedure is used, the recombination procedure 
and separation cuts influence the resulting unweighted events. 

PRECOMB (1:4): recombination parameters: 

PRECOMB (1) : separation cut on angle between beam and photon in degrees if 
SRECOMB = 1,2,3, 

PRECOMB (2) : recombination cut for photon energy in GeV if SRECOMB = 1, 2, 3, 

PRECOMB (3) : invariant-mass recombination cut in GeV if SRECOMB = 1,3 or angular 
recombination cut for leptons in degrees if SRECOMB = 2, 

PRECOMB (4) : angular recombination cut for quarks in degrees if SRECOMB = 2. 
The recombination scheme is defined in the subroutine RECOMBINATION in public. f. 
The recombination scheme for SRECOMB = 1 is described in detail in Ref. [6] for the 
parameters PRECOMB(l) = 5, PREC0MB(2) = 1, and PREC0MB(3) = 5 or 25, called "bare" 
or "calo" cuts, respectively. For all three photon recombination schemes (SRECOMB = 
1,2,3), the photon momentum is set to zero if the angle between the photon and the 
beam is smaller than PRECOMB (1). This cut is applied at the very beginning before 
photon recombination is performed. 

For SRECOMB = 1, the minimal invariant mass of the photon and any charged final- 
state fermion, M(7,FC0MB), is determined, and the photon is recombined with the 
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final-state fermion FCOMB if M(j, FCOMB) < PREC0MB(3) or the photon energy < 
PREC0MB(2). Explicitly, this means that the photon momentum is added to the momen- 
tum of the fermion FCOMB, and the photon is discarded. 

For SRECOMB = 2, the minimal angle between the photon and any charged final-state 
fermion, #(7, FCOMB), is determined, and the photon is recombined with the fermion 
FCOMB if 6'(7,/ ± ) < PREC0MB(3) in case / ± is charged lepton of the final state or if 
$(7? <?) < PREC0MB(4) in case q is a final-state quark, or if the photon energy E 1 < 
PREC0MB(2). 

For SRECOMB = 3 the minimal invariant mass of the photon and any charged fermion 
(including the incoming electron and positron), M("f, FCOMB), is determined, and the 
photon is recombined with the fermion FCOMB if M(7, /) < PREC0MB(3), where / is one 
of the initial- or final-state charged fermions, or if E 1 < PREC0MB(2). In case FCOMB is a 
beam electron or positron, the photon momentum is set to zero. 

SQCDEPEM: flag for the inclusion of gluon-exchange diagrams in the tree-level mode: 
SQCDEPEM = : gluon-exchange background diagrams are not included, 
SQCDEPEM = 1 : gluon-exchange background diagrams are included, 
SQCDEPEM = 2 : only gluon-exchange background diagrams are included. 
The choice of SQCDEPEM only affects the tree-level cross sections to e + e~ — > 4/ and 
e + e~ — > 4/7, and for SISR = 1 also the LL corrections from ISR. This QCD contribution 
then consists of gluon-exchange background diagrams to the tree-level processes e + e" — > 
4/ and e + e _ — > 4/7 as described in Ref. [4]. In the RC mode gluon-exchange diagrams 
can be taken into account only in the tree- level matrix element to e + e~ — + 4f. 

SATGC: flag for the inclusion of anomalous triple gauge-boson couplings (ATGCs): 
SATGC = : ATGCs are not included, 
SATGC = 1 : ATGCs are included. 

If SATGC = 1, the parameters GTGC(1 : 11 , 1 : 2) must be specified in the input file. They 
are stored in the common block ATGC and correspond to the ATGCs AgY , Ak, v , X v , 
k v , ~\ v , gX$, f¥ >5 , h\ fi with V = 7(= A), Z. The ATGCs affect only the tree-level cross 
section for e + e~ — > 4/, and for SISR = 1 also the LL corrections from ISR are applied 
to the tree-level cross sections involving ATGCs. 

QAQGC: flag for the inclusion of anomalous quartic gauge-boson couplings (AQGCs): 
QAQGC = : AQGCs are not included, 
QAQGC = 1 : AQGCs are included. 

If QAQGC = 1, the parameters ALAM2(0:4) must be specified in the input file. They are 
stored in the common block AQGC and correspond to the AQGCs ao iC ,n/A 2 and ao, n / 'A 2 . 
The inclusion of AQGCs is not supported in the RC mode and only affects the results 
obtained in the tree-level mode (SRC = 0, SB0RN5 = 1) and the IBA mode (SRC = 3) 
for e + e~ — ■> 4/7. 

SCUTS : flag for choosing the separation cuts: 
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SCUTS = : no separation cuts are applied, 

SCUTS = 1 : ADLO cuts are applied as defined in racoonww.f and in Ref. [10], 

SCUTS = 2 : LC cuts are applied as defined in racoonww.f, 

SCUTS = 10 : a minimal set of cuts must be specified in the input file, 

SCUTS = 11 : cuts on the energies of all final-state particles, on all angles, and on 

invariant masses must be specified in the input file. 
All energies and invariant masses are measured in GeV and all angles in degrees. Non- 
standard cuts, i.e. cuts which are not covered by the choices of SCUTS that RacoonWW 
provides, can be included at the end of the subroutine CUT in public.f. 

4-2 Options not provided in the input file 

RacoonWW has a few more options that are not provided in the input file. Changing 
those options is only advisable if the user is familiar with the details of the calculation 
as described in Ref. [6]. If they are common to both the slicing and subtraction branches, 
RACOONWW sets them to their default values in subroutine INITIALIZE in racoonww.f 
(after "setting options"). Special options for the slicing and subtraction branches are 
set in the subroutines INITSLICING and INITSUBTRACTION, respectively. 

The following flags are set in subroutine INITIALIZE in racoonww.f: 

SON : flag for the projection of the phase space of the final-state particles to the one with 
on-shell W bosons. RacoonWW provides four different versions for the projection of 
off-shell momenta to on-shell momenta: 

SON = : directions of momenta of W + , / 3 , and / 5 are fixed (default), 

SON = 1 : directions of momenta of W + , / 4 , and / 5 are fixed, 

SON = 2 : directions of momenta of W + , / 3 , and / 6 are fixed, 

SON = 3 : directions of momenta of W + , and f§ are fixed 

for the process e + e~ -> f^Uhh with W+ ~> hh and W ~ ~> hh- 
Since the projections are performed in the CM frame of the incoming e + e~, the di- 
rection of the W~ is opposite to the direction of the W + . The projection is done in 
the subroutine 0FT00N in eeWW4f_DPA.f . In Ref. [6] the accuracy of the DPA has been 
investigated by modifying the implementation of the DPA and comparing the results. 
The comparison of the on-shell projections, SON = 1 ("def") and SON = 2 ("proj"), is 
one of the considered options. 

QFAST : flag for the evaluation of the virtual corrections in DPA: 

QFAST = : form factors for virtual corrections from full formulas (slow), 

QFAST = 1 : evaluation from expansion into Legendre polynomials (fast); this is the 

default value of QFAST. 
Details of this approach can be found in Section 3.2.3 and in Ref. [6]. 

QSOFT: flag for the definition of the DPA part of radiative corrections. Terms that are 
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subtracted from the virtual corrections in DPA are fixed by 

QSOFT = 1 : endpoint part of Ref. [19] as given in Eq. (4.29) of Ref. [6] (default), 

QSOFT = 2 : only logarithms as defined in Eq. (4.2.53) of Ref. [20], 

QSOFT = 3 : YFS factor for virtual photons as defined in Eq. (4.31) of Ref. [6]. 
Note that results obtained with different choices for QSOFT differ only by terms that are 
beyond DPA accuracy. In Ref. [6] the accuracy of the DPA has been investigated by 
modifying the implementation of the DPA and comparing the results. The definition of 
the finite virtual corrections, QSOFT = 1 ("def") and QSOFT = 3 ("eik"), is one of the 
considered options. 

RacoonWW provides the possibility to use PAW/HBOOK [45,46] for histogram filling 
and plotting. The usage of PAW/HBOOK is only allowed when generating weighted 
events (NEVENTSUNW=0). This option can be activated as follows: 

• Set SPAW = 1 in the subroutine INITIALIZE in racoonww.f. 

• Uncomment all lines in the if-environment starting with "if (spaw. eq. 1. . ." in the 
subroutine INITIALIZE. 

• Edit the makefile and add pawgraphs.o to the list of the objects OBJS. When using 
HBOOK/PAW, CERNLIB must be linked in the makefile. 

In pawgraphs . f the histograms are defined in the subroutine SETUP -GRAPHS and filled in 
the subroutine GRAPHS, which are called in kern.f. This is the place where the user can 
add and modify histograms. 

In the following we describe some branch-specific options of RacoonWW. 



4-2.1 Options specific for the slicing branch 

All the flags and parameters specific for the slicing branch are set in the subroutine 
INITSLICING in slicing. f: 

OPT : flag for the optimization of the a priori weights: 
OPT = : no optimization, 

OPT = 1 : self-optimization of the a priori weights (default). 

In the slicing branch the self-optimization of the a priori weights is performed when 
OPT = 1 in slicing.f. The number of events used for the optimization of the a priori 
weights and how often the optimization is performed can be separately specified for the 
calculation of the 4/ weights (N0PT,I0PT) and the 4/7 weights (NOPTG, I0PTG). The 
default setting is NOPT = 10 6 , I0PT = 10, and NOPTG = 25 x 10 6 , I0PTG = 1. The flag 
OPT and the parameters enter in function KERN in kern.f. The a priori weights are 
optimized as long as the number of generated events is smaller than NOPT or NOPTG. For 
the e + e~ — > 4/ part the a priori weights are optimized after 5 x 10 3 , 10 x 10 3 , 15 x 10 3 , 
25 x 10 3 generated events and when the number of events, NEVENTS, is a multiple of 
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NOPT/IOPT as long as NEVENTS < NOPT. For the e + e~ — »■ 4/7 part the a priori weights 
are optimized after 5 x 10 5 , 15 x 10 5 , 30 x 10 5 generated events and when NEVENTS is 
a multiple of N0PTG/I0PTG as long as NEVENTS < NOPTG. 



4-2.2 Options specific for the subtraction branch 

The flags and parameters specific for the subtraction branch are set in the subroutine 
INITSUBTRACTION in subtraction . f: 

SSUB : flag for the definition of the subtraction functions: 
SSUB = 1 : subtraction terms as defined in Ref. [19], 
SSUB = 2 : subtraction terms as defined in Ref. [20]. 

The results obtained by the two subtraction variants differ only by terms that are 
beyond the DPA accuracy. More details can be found in Ref. [6]. 

NOPT : parameters for the adaptive optimization. 

NOPT defines the numbers of events after which the adaptive optimization is performed, 
i.e. after which event the a priori weights are recalculated. The optimization steps are 
given by N0PT(I, NG) x number of channels, where NG = 1, 2, 3 denotes the phase-space 
generators discussed in Section 3.1.5. The adaptive optimization stops after the i-th 
optimization step if N0PT(I + 1, NG) = 0. The default settings are N0PT(I, NG) = I x 100 
for I = 1, . . . , 8 and zero otherwise. 

ALPHAMIN: the a priori weights have a minimal value which is fixed by setting the pa- 
rameter ALPHAMIN. 



4.3 The output of Racoon WW 



The output consists of an output file and 27 (or more, if the user added histograms) 
data files, which contain the data for the distributions. The name of the output file is 
set in the input file and stored in the variable 0UTPUTFILE. For SPAW = 1 there is also a 
PAW file pawplot.paw with the histograms generated with PAW/HBOOK [45,46]. If 
I0UT 7^ is set in the subroutine INITIALIZE in racoonww.f, a file optimization, info 
will be generated. This file contains information about the kinematical channels used 
in the multi-channel integration and about the optimization of the a priori weights as 
described in Section 3.1. By default I OUT = is set. 

The event information is stored in the common block EVENT (see Section 2). The weights 
of the slicing and subtraction branches are described in Table 5 for SMC = 1, in Table 6 
for SMC = 2, and in Table 7 for SMC = 3. 
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WEIGHT (1) : tree-level cross section to e + e" -> 4/, vanishes for SBDRN4 = or SISR = 1, 
WEIGHT (2) : e+e" -> WW -> 4/ DPA part of 0(a) corrections, 

WEIGHT (3) : virtual+soft+fmal-state collinear photon part (to avoid double counting, LL 
0(a) corrections from ISR are subtracted when SISR = 1), 

WEIGHT (4) : initial-state collinear photon part of 0(a) corrections (to avoid double count- 
ing, LL 0(a) corrections from ISR are subtracted when SISR = 1), collinear 
photon radiation off e + , 

WEIGHT (5) : same as WEIGHT (4) but for collinear photon radiation off e", 

WEIGHT (6) : SRC = 1 : the tree-level cross section is convoluted with the ISR structure func- 
tions including LL corrections up to 0(a 3 ) 

SRC = 2 : IBA for e+e" -► WW -► 4/, the CC03 IBA cross section is convo- 
luted with the ISR structure functions, 

vanishes for SISR = 0, 

WEIGHT (7) : SRC = : tree- level process e+e~ -> 4/7, 

SRC = 1 : e + e~ — ► 4/7 bremsstrahlung contribution, 

SRC = 3 : IBA for e+e" -► 4/7. 



Table 5 

Description of the weights WEIGHT (1 :7) of the slicing branch (SMC = 1) 



WEIGHT (1) 


: tree- level cross section to e + e~ — ► 4f, 




vanishes for SBDRN4 = or SISR = 1, 


WEIGHT (2) 


: SRC = 0, 1 : tree-level cross section for e + e~ — ► 4/, 




SRC = 2 : IBA cross section for e+e~ -> WW -> 4/, 




convoluted with ISR structure functions, LL corrections from ISR up to 




0(a 3 ), collinear photon radiation off e + and e~, 




vanishes for SISR = 0, 


WEIGHT (3: 5) 


: 4/ part of radiative corrections including virtual 0(a) corrections and 




subtraction functions for SRC = 1 (0(a) LL corrections from ISR omitted 




to avoid double counting), 




vanishes for SRC = 0, 2, 3, 




WEIGHT (3) : contribution with 4/ kinematics, 




WEIGHT (4) : collinear photon radiation off e + , 




WEIGHT (5) : collinear photon radiation off e~, 


WEIGHT (6) 


: SRC = : tree- level process e + e~ — ► 4/7 when SB0RN5 = 1, 




SRC = 1 : bremsstrahlung part of 0(a) corrections, 




SRC = 3 : IBA for e+e" -► 4/7, 


WEIGHT (7) 


: presently not used, 


WEIGHT (8: 29) 


: contributions from the subtraction functions to the 4/7 part for SRC = 1, 




vanishes for SRC = 0,2,3. 



Table 6 

Description of the weights WEIGHT (1 :29) of the subtraction branch (SMC = 2) 
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WEIGHT (1) : tree-level cross section to e+e" -> 4/, e+e" -> WW -> 4/ DPA part of 0(a) 
corrections, virtual+soft+final-state collinear photon part, and soft part of 
the 0(a 2 ) LL corrections from ISR when SISR = 1, 

WEIGHT(2:3) : vanishes, 

WEIGHT (4) : initial-state collinear photon part of 0(a) corrections, collinear part of the 
0(a 2 ) LL corrections from ISR when SISR = 1, collinear photon radiation 
off e+, 

WEIGHT (5) : same as WEIGHT (4) but for collinear photon radiation off e", 

WEIGHT (6) : collinear part of the 0(a 2 ) LL corrections from ISR when SISR = 1, collinear 

photon radiation off both e + and e~ simultaneously, 
WEIGHT (7) : e + e~ — > 4/7 bremsstrahlung contribution. 



Table 7 

Description of the weights WEIGHT(1:7) of the slicing branch (SMC 
supported in the RC mode (SRC = 1) 



3). This branch is only 



4.3.1 The output file of Racoon WW 



In the output file RacoonWW provides the following information: 

• the input parameters, the effective branching ratios, and the choices for the options, 

• the total cross section together with the maximal values of the weights (in the subtrac- 
tion branch also the corresponding event numbers and channels are given), and 

• the total cross sections for the subcontributions (tree-level 4/ cross section, 2^4 part 
of the corrections, ISR up to 0(a 3 ), and 2-^5 part of the corrections) together with 
the maximal values of the weights. 

Finally, the number of rejected events is given. Events are discarded owing to the applied 
separation cuts and in case of the slicing branch also owing to the imposed technical cuts. 
RacoonWW provides the number of rejected events in the calculation of each weight 
separately. 



4-3.2 The data files for distributions 

The data for the distributions are written in files dat . 01,. . .,dat . 27. The histograms 
are defined in the subroutine SETTINGS in public. f as described in Table 8. The file 
dat. 01 contains the total cross section for test purposes. Note that the momenta of the 
initial-state particles after ISR, P(l : 2, J,K), do not correspond to physically measurable 
quantities and, thus, must not be used for defining histograms. 

The data files contain the observable (column 1), the differential cross section for the 
chosen setup (column 2) and the corresponding lowest-order cross section of the tree-level 
process e + e~ — > 4/ (column 4). The statistical errors are provided in columns 3 and 5, 
respectively. These are calculated as described in Eq. (15) and below. The histograms are 
normalized such that the total cross section is obtained as a = Sbins bin-height x bin-width. 
The angles of the histograms dat. 21 and dat. 22 are denned in Eqs. (7.4) and (7.3) of 
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data file 


distribution in the 


dat 


02 


invariant mass of fa and /4, [75:85] GeV (in case of CC processes this corresponds 






to the W + invariant mass), 




dat 


03 


invariant mass of f$ and f$, [75:85] GeV (in case of CC processes this corresponds 






to the W~ invariant mass), 




dat 


04 


cosine of the angle between V\ and /i = e + , [— 1:1], 




dat 


05 


cosine of the angle between V" 2 and ft = e~, [— 1:1], 




dat 


06 


cosine of the angle between V\ and f%, [—1:1], 




dat 


07 


cosine of the angle between V\ and /4, [—1:1], 




dat 


08 


cosine of the angle between Vi and /s, [—1:1] 




dat 


09 


cosine of the angle between V2 and f$, [—1:1], 




dat 


10 


energy of the photon, [0:50] GeV (4/ weights are added to the first bin), 




dat 


11 


cosine of the angle between 7 and j\ = e + , [—1:1] (4/ weights are added to the 






last bin), 




dat 


12 


minimal angle between 7 and charged final-state fermions, [0:180] degrees 
weights are added to the first bin), 


(4/ 


dat 


13 


energy of / 3 , [0:^/2] GeV, 




dat 


14 


energy of / 4 , [0:^/2] GeV, 




dat 


15 


energy of / 5 , [0:^/2] GeV, 




dat 


16 


energy of / 6 , [0:^/2] GeV, 




dat 


17 


cosine of the angle between V\ and in the rest frame of V\, [—1:1], 




dat 


18 


cosine of the angle between V\ and in the rest frame of Vi, [—1:1], 




dat 


19 


cosine of the angle between V2 and /s in the rest frame of V2, [—1:1], 




dat 


20 


cosine of the angle between V2 and Jq in the rest frame of V2, [—1:1], 




dat 


21 


angle between decay planes of ^3, J4 and /s, f$, [0:360] degrees, 




dat 


22 


angle between decay planes of V\, f\ = e + and V±, /4, [0:360] degrees, 




dat 


23 


angle between decay planes of V2, /b = e~ and V2, /s, [0:360] degrees, 




dat 


24 


cosine of the angle between ^3 and f^, [—1:1], 




dat 


25 


cosine of the angle between f$ and f$, [—1:1], 




dat 


26 


cosine of the angle between V\ and j\ = e + in the CM frame of V\ and V2, [— 


1:1], 


dat 


27 


cosine of the angle between V2 and / 2 = e~ in the CM frame of V\ and V2, [- 


1:1]. 



Table 8 

Description of the data files for the differential cross sections to /1/2 — ► V1V2 — > fsfifefeil) 
where f\ = e + , / 2 = e~, V\ — > fzh, and V2 — > /s/6 and the e + beam goes in the positive 
z direction. All histograms are defined with 50 bins. 
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Ref. [6], respectively. 

In the subroutine SETTINGS in public. f the user can define additional histograms as 
described at the end of SETTINGS. The histograms are filled in 

SUBROUTINE HISTOGRAM (MIN , MAX , X , Y , HISTO , N , NAME , PARTS , STEPS) 

MIN,MAX : minimal and maximal value of the X range, 

X,Y : X and Y value, Y is the weight of the event, 

HISTO : number of the histogram (HISTO < MAXH), 

N : total number of (weighted) events, used only for STEPS = 4, 

NAME : name of the data file, 

PARTS : number of bins of the histogram (PARTS < MAXP) , 
STEPS : initialization, filling, and output of the histogram, 
STEPS = 1 : initialization of the histogram, 

STEPS = 2 : filling the histogram with one of the weights of an event (column 
2), " 

STEPS = 3 : evaluation of the histogram information for a complete event 
(column 2,3), and filling the histogram with a single event of 
the tree-level cross section to e + e~ — > 4/ (column 4,5), 

STEPS = 4 : writing data file of the histogram 

at the end of public. f. The parameters MAXH and MAXP specify the maximal number of 
histograms and bins and are set to MAXH = 50 and MAXP = 1000, respectively. If more 
histograms or bins are required, these parameters have to be increased. 

For the subtraction branch the splitting into STEPS = 2 and STEPS = 3 is needed since the 
weights for the bremsstrahlung process e + e~ — > 4/7 and the corresponding subtraction 
functions are related to different phase-space points and, hence, cannot be combined into a 
single weight. In the soft and collinear limits the momenta corresponding to these different 
weights approach each other and the weights approximatively cancel, leading to IR-safe 
observables. These cancellations of the soft and collinear singularities occur within single 
histogram bins. Therefore, the histograms must be filled in a first step (STEPS = 2) with all 
subcontributions corresponding to one weighted event. From the resulting contributions 
to each histogram bin the corresponding statistical errors can be calculated (STEPS = 3). 

In addition, for SPAW = 1 the HBOOK file pawplot .paw is generated. It contains the his- 
tograms for the distributions as specified in pawgraphs.f. In racoonww.f the subroutine 
SETUP _GRAPHS is called where the histograms are defined (HBOOK). The subroutine GRAPHS 
is called in racoonww.f after the weights have been calculated. There the observables 
are defined and the corresponding histograms are filled (HFILL). In addition three more 
subroutines, GRAPHS_B, GRAPHS_K, and GRAPHS_R, are provided which generate histograms 
showing the distributions separately for the tree-level e + e~ — >■ 4/ contribution, the whole 
4/ part, and the 4/7 part. As usual, the histograms can be viewed by opening the PAW 
[45] main browser, and loading the file pawplot. paw. 
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4-4 The unweighting procedure and the interface with Pythia 

As default, RACOON WW generates only weighted events. For SRC = 1 and SMC = 3, or for 
SRC = 0, 2, 3, RacoonWW can generate unweighted events. The unweighting procedure 
is activated as follows: 

• The number of generated unweighted events NEVENTSUNW must be set in the main pro- 
gram in racoonww.f . By default NEVENTSUNW = 0, and no unweighted events are gener- 
ated. In the RC mode, to obtain a precision comparable to the "best" result of Ref. [6], 
we recommend to use NEVENTSUNW = 5 x 10 4 . With an efficiency of one unweighted event 
per 10 3 weighted events, this corresponds to the calculation of approximately 5 x 10 7 
weighted events. In the tree-level and IBA modes the efficiency is about a factor 100 
better and we recommend to use NEVENTSUNW = 5 x 10 5 . Note that for certain processes 
the efficiency can strongly depend on the separation cuts. 

• The number of weighted events (NEVENTSW) that are used only to determine the maximal 
weights for the hit-and-miss algorithm is set in the input file. This should be much 
smaller than what is usually used without unweighting, for instance, we recommend 
NEVENTSW = 10 5 . 

When including radiative corrections (SRC = 1), only the branch SMC = 3 can be used in 
the unweighting procedure. 

The unweighted events are generated in the subroutine EVENTGENERATION in racoonww.f 
from weighted events with a hit-and-miss algorithm. More precisely, for each weighted 
event an unweighted event is generated if the corresponding weight is larger than the 
product of a random number RAND (0 < RAND < 1) and the maximal weight. Other- 
wise no unweighted event is generated. The maximal weight, WEIGHTT0TMAX, is calcu- 
lated by multiplying the maximal weight of the first NEVENTSW events with a factor two. 
The number NEVENTSW and the factor can be tuned to suppress weights that are larger 
than WEIGHTT0TMAX in the main calculation. If a weight exceeds the maximal weight, 
WEIGHTT0TMAX is redefined by the value of the weight and a warning is written in the 
output file. 

The maximal weight of the bremsstrahlung contribution exceeds by far the maximal 
weights of the other subcontributions. Therefore the other weights are not calculated 
for each event in order to enhance the efficiency of the unweighting procedure. They 
are calculated only all NWEIGHT times but are multiplied by NWEIGHT. In this way the 
probability of an unweighted event from the different subcontributions remains unchanged. 
The integers NWEIGHT (1 : 29) are chosen in such a way that each subcontribution leads to 
roughly the same maximal weight, which has to be smaller than WEIGHTT0TMAX. 

The information about the (unweighted) event is provided in the HEP standard format 
in the common block HEPEVT [24] . The parameters and flags of the common block HEPEVT 
are set in the subroutine EVENTGENERATION in racoonww.f. For a detailed description we 
refer to Ref. [23], p. 63. While the first six particle entries in HEPEVT are the external 
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fermions, e + e~ — > 73/4/5/6, the last up to two entries are filled by a possible hard photon 
from the bremsstrahlung process e + e~ — > 4/7 and, if higher-order LL corrections from 
ISR are included, a photon collinear to the beam. The four-momentum of the ISR photon 
collinear to the beam is calculated from the four-momenta of the other external particles 
using momentum conservation. Note that the positron beam goes in +z direction and 
that cuts and photon recombination influence the resulting unweighted events. 

Although negative weights appear in the RC mode also for SRC = 3, they have only a negli- 
gible influence on the observables that we have studied so far. Nevertheless, RacoonWW 
also unweights the events with negative weights, but does not fill the HEPEVENT common 
block with these weights. The resulting unweighted negative- weight events can, for in- 
stance, be used to estimate the uncertainty resulting from discarding negative weights. 

In the output file, RacoonWW provides the number of unweighted events, the number 
of unweighted events with negative weights, and the number of events that enforce a 
redefinition of the maximal weight. Note that the result for the total cross section and for 
the distributions are always calculated from the weighted events. 

To include parton shower and hadronization for coloured final states, RACOON WW pro- 
vides an interface with Pythia for unweighted events by calling the Pythia subroutine 
PY4FRM (see p. 236 of Ref. [23]): 

SUBROUTINE PY4FRM ( ATOTSQR , A1SQR , A2SQR , ISTRAT , IRAD , ITAU , ICOM) 

ATOTSQR : total matrix-element squared for the event calculated in the subroutine 
EVENTGENERATION as ATOTSQR = ATOTSQ(NS), where ATOTSQ(NS) is the ma- 
trix element squared of the NS-th weight WEIGHT (NS), 

A1SQR : matrix-element squared for the configuration when fa, and /s, fe are colour 
singlets (A1SQR = AISQ(NS)), 

A2SQR : matrix-element squared for the configuration when / 3 , / 6 and /s, f'4 are colour 
singlets (A2SQR = A2SQ(NS)), 

ISTRAT : flag for the strategy to select either of both colour configurations; we choose 
ISTRAT = 0, 

IRAD : flag for final-state photon radiation in Pythia. As default, photon radiation 
is switched off (IRAD = 0) to avoid double counting in the RC mode (SRC = 
1); in the tree- level and IBA modes IRAD = 1 may be chosen, 

ITAU : flag for the handling of r-lepton decays; by default r decays are not included 
(ITAU = 0), since the final-state fermions are considered to be massless in 
RacoonWW; however, RacoonWW provides a subroutine to generate a 
massive two-particle phase space (see below), so that r decays can be taken 
into account (ITAU =1), 

ICOM : flag for the place where the information about the event is stored; by default 
ICOM = in order to store the information in the HEPEVT common block. 

The amplitudes AT0TSQ(1:8), A1SQ(1:8), A1SQ(1:8), corresponding to the weights in 
Tables 5, 6, and 7, are calculated by using AT0TSQ4, A1SQ4, A2SQ4 for the tree-level pro- 
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cesses e+e" -> 4/ and by AT0TSQ5, A1SQ5, A2SQ5 for the tree-level processes e+e" -> 4/7. 
The colour amplitudes are calculated in the functions M2_EPEM_4F and M2_EPEM_4FA in 
ee4fa_amps.f and stored in the common blocks COLAMP and COLAMPA. 

As default, this interface is not active in RacoonWW. It can be activated as follows: 

• Link Pythia by adding the file with the Pythia code, e.g. pythial568.f of version 
6.1, to the list of objects OBJS in the makefile. Pythia can be downloaded from 
http:/ /www. thep.lu.se/| ~ torbjorn/Pythia.html. 

• Uncomment the line with the call of PY4FRM in the main program in racoonww.f . 

Note that the compilation of Pythia needs extensive memory. 

As mentioned above, RacoonWW provides the possibility of converting pairs of massless 
momenta to massive momenta by calling 

SUBROUTINE MASSIVE (Ml , M2 , PHEP1 , PHEP2) 

M1,M2 : masses of the two final-state particles, 

PHEP1,PHEP2 : four momenta of the two final-state particles with zero masses before 
the call of MASSIVE and with masses Ml ,M2 after the call of MASSIVE 

in racoonww . f . This could be used to switch to massive kinematics for the r lepton before 
interfacing with Pythia. Note that the momenta of two final-state particles have to be 
changed in order to maintain four-momentum conservation and the mass-shell conditions. 



5 Summary 



We have described the MC generator RacoonWW, version 1.3, for four-fermion produc- 
tion, off-shell W-pair production, as well as the radiative processes e + e~ — > 4/7. At tree 
level, the complete set of diagrams is included. For e + e~ — * WW — * 4/, the prediction is 
based on the complete 0(a) corrections in double-pole approximation and QED initial- 
state radiation up to 0(a 3 ). The corresponding theoretical uncertainty is at the level of 
0.5% for the total cross section and of the order of 1% for distributions in the energy 
range between ~ 170 and 500 GeV. The use of a double-pole approximation, in general, 
puts a restriction on the theoretical uncertainty at the level of a few per mil away from 
threshold. RacoonWW also provides an improved-Born approximation with an uncer- 
tainty of 2%, which is also applicable in the threshold region. More precise calculations 
require the inclusion of radiative corrections beyond double-pole approximation. 
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A Sample runs and output 

Two test runs are discussed in the following. One is a run of RacoonWW in the RC 
mode. The second one is the calculation of the cross section to e + e~ — > 4/7 within the IBA 
mode. Both test runs and three additional sample runs are discussed in the file README. 

A.l Test run in the RC mode 

The first test run calculates the "best" predictions for the cross section to the process 
e + e~ — > ud/i~z/ M (7) for a CM energy of s/s = 200 GeV. Thereby, the full 0(a) corrections 
in DPA and leading higher-order corrections are included, and the W-boson width is 
calculated. The recombination scheme and separation cuts of Ref. [6] are applied. The 
numerical integration is performed in the slicing branch with 5 x 10 7 (weighted) events. 
On a DEC/ ALPHA workstation the computing time is about 30 hours. 

A. 1.1 Input file input sli 



outputfile ! name of output file 

200d0 ! energy: CMF energy (in GeV) 

50000000 ! neventsw: number of weighted events 

1 ! smc: choice of MC branch: l(or 3): slicing 2 : subtraction 

1 ! sborn4: include Born ee->4f : 0:no 1-3: yes 

1 ! sborn5: include Born ee->4f +photon: 0:no l:yes 

! sborng5: include Born ee->4f +gluon: 0:no l:yes 

1 ! sisr: include higher-order ISR: 0:no l:yes 

1 ! src: include radiative corrections: 0:no 1:DPA 2:IBA-4f 3:IBA-4fa 
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: 1 VaC U-LXLlc! XdillULLCL ri 


OHO 


1 Tfir 1-i Trip 1amh<-ia~7 


OHO 


! TGC f 4~A 


flHO 


! 1 \J\j X Zj 


UaU 


! IbC I_D A 


OdO 


! TGC f_5~Z 


OdO 


! TGC h_l~A 


OdO 


! TGC h_l~Z 


OdO 


! TGC h_3"A 


OdO 


! TGC h_3~Z 
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! qaqgc: anomalous quartic gauge couplings (QGC) : 0:no l:yes 

OdO ! QGC a_0/Lambda~2 

OdO ! QGC a_c/Lambda~2 

OdO ! QGC a_n/Lambda~2 

OdO ! QGC tilde a_0/Lambda~2 

OdO ! QGC tilde a_n/Lambda~2 

10 ! scuts: separation cuts: 0:no l,2:default(ADL0,LC) 10,ll:input 

OdO ! photon (gluon) energy cut 

OdO ! charged-lepton energy cut 

OdO ! quark energy cut 

OdO ! quark-quark invariant mass cut 

OdO ! angular cut between photon and beam 

OdO ! angular cut between photon and charged lepton 

OdO ! angular cut between photon (gluon) and quark 

OdO ! angular cut between charged leptons 

OdO ! angular cut between quarks 

OdO ! angular cut between charged lepton and quark 

lOdO ! angular cut between charged lepton and beam 

lOdO ! angular cut between quark and beam 



A. 1.2 Output file 

smc= 1 : Phase-space-slicing branch of RacoonWW 



technical cutoff parameters (photon) : 
delta. s = 1.000000000000000E-003 
delta.c = 5.000000000000000E-004 

Input parameters: 

CMF energy = 200.00000 GeV, Number of events = 50000000, 

alpha(O) = 1/ 137.0359895, alpha(MZ) = 1/128.88700, alpha. s = 0.11900, 

GF = .1166370E-04, 

MW = 80.35000, MZ = 91.18670, MH = 150.00000, 

GW = 2.08699, GZ = 2.49471, 

me = .51099907E-03, mmu = 0.105658389, mtau = 1.77705, 

mu = 0.00485, mc = 1.55000, mt = 174.17000, 

md = 0.00485, ms = 0.15000, mb = 4.50000. 

Effective branching ratios: 

leptonic BR = 0.32520, hadronic BR = 0.67480, total BR = 1.00000 
Process: anti-e e -> u anti-d mu anti-nu_mu (+ photon) 
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pp- 


o 

V . 


pm= 


0. 


qalp= 


2: 


qgw= 


4: 


nnrnn = 


J. ■ 


sborn4= 


1: 


ssigepem4= 


0: 


qqcd= 


2: 


src= 


1: 


ssigepem5= 


n . 

t> . 


qqcd= 


2: 


qreal= 


1: 


qnf = 


3: 


sisr= 


1: 


srecomb= 


1: 



one-loop W-boson width calculated (with QCD corr.) 



tree-level process ee -> 4f . 

all electroweak diagrams included. 

naive QCD corrections included. 

virtual corrections in DPA and real corrections included, 
real photon corr. : only CC11 class of diagrams included, 
naive QCD corrections included 

imaginary part of virtual corrections neglected. 

off -shell Coulomb singularity with off -shell Born included. 

initial-state radiation up to order alpha~3 included. 



precomb(l) = 

precomb(2) = 

precomb(3) = 

precomb(4) = 



5.00000 
1.00000 
5.00000 
0.00000 



scuts=10 : 

angle (1,3) > 

angle (1,4) > 

angle (1,5) > 

angle (2, 3) > 

angle (2,4) > 

angle (2, 5) > 



with separation cuts: 
10.00000 deg 
10.00000 deg 
10.00000 deg 
10.00000 deg 
10.00000 deg 
10.00000 deg 



events 


intermediary 


results 


preliminary 


results 


1000000 


592 


61877 


+- 


2 


47251 


592 


61877 


+- 


2 


47251 


2000000 


591 


72383 


+- 


1 


89683 


592 


17130 


+- 


1 


55815 


3000000 


594 


65445 


+- 


1 


89273 


592 


99902 


+- 


1 


21535 


4000000 


593 


26888 


+- 


1 


80112 


593 


06648 


+- 


1 


01667 


5000000 


587 


37590 


+- 


1 


65706 


591 


92837 


+- 





87826 


6000000 


591 


10834 


+- 


2 


22044 


591 


79170 


+- 





82013 


7000000 


594 


42877 


+- 


1 


76546 


592 


16842 


+- 





74684 


8000000 


590 


00165 


+- 


1 


68242 


591 


89757 


+- 





68649 


9000000 


589 


54677 


+- 


1 


59058 


591 


63637 


+- 





63529 


10000000 


594 


99013 


+- 


1 


90289 


591 


97175 


+- 





60260 


11000000 


588 


48466 


+- 


1 


82284 


591 


65474 


+- 





57233 


12000000 


593 


04183 


+- 


1 


90019 


591 


77033 


+- 





54801 


13000000 


592 


36715 


+- 


1 


68580 


591 


81624 


+- 





52221 
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14000000 


590 


27937 


+- 


1 


64218 


591 


70646 


+- 





49890 


15000000 


591 


71289 


+- 


1 


79049 


591 


70689 


+- 





48070 


16000000 


591 


53899 


+- 


1 


76907 


591 


69640 


+- 





46402 


17000000 


587 


50308 


+- 


1 


58333 


591 


44973 


+- 





44654 


18000000 


587 


89929 


+- 


1 


55550 


591 


25249 


+- 





43050 


19000000 


592 


93474 


+- 


1 


73210 


591 


34103 


+- 





41790 


20000000 


592 


16324 


+- 


1 


74407 


591 


38214 


+- 





40647 


21000000 


592 


24008 


+- 


1 


80517 


591 


42299 


+- 





39655 


22000000 


592 


44203 


+- 


1 


72115 


591 


46931 


+- 





38652 


23000000 


588 


26438 


+- 


1 


68691 


591 


32997 


+- 





37692 


24000000 


590 


49483 


+- 


1 


81216 


591 


29517 


+- 





36902 


25000000 


591 


15275 


+- 


1 


79673 


591 


28947 


+- 





36148 


26000000 


593 


47325 


+- 


1 


78914 


591 


37346 


+- 





35432 


27000000 


591 


22938 


+- 


1 


63230 


591 


36813 


+- 





34651 


28000000 


591 


06261 


+- 


1 


77909 


591 


35722 


+- 





34013 


29000000 


590 


90543 


+- 


1 


71644 


591 


34164 


+- 





33369 


30000000 


591 


59423 


+- 


1 


68095 


591 


35006 


+- 





32740 


31000000 


591 


88912 


+- 


1 


70980 


591 


36745 


+- 





32160 


32000000 


588 


95475 


+- 


1 


67046 


591 


29205 


+- 





31589 


33000000 


589 


07952 


+- 


1 


59758 


591 


22500 


+- 





31012 


34000000 


591 


09482 


+- 


1 


83148 


591 


22117 


+- 





30578 


35000000 


593 


56319 


+- 


1 


67756 


591 


28809 


+- 





30089 


36000000 


592 


41598 


+- 


1 


67903 


591 


31942 


+- 





29623 


37000000 


593 


56115 


+- 


1 


88861 


591 


38001 


+- 





29270 


38000000 


588 


14474 


+- 


1 


69946 


591 


29487 


+- 





28849 


39000000 


591 


59522 


+- 


1 


65363 


591 


30257 


+- 





28427 


40000000 


590 


01906 


+- 


1 


63315 


591 


27048 


+- 





28016 


41000000 


590 


88361 


+- 


1 


92445 


591 


26105 


+- 





27732 


42000000 


592 


07152 


+- 


1 


69336 


591 


28034 


+- 





27371 


43000000 


590 


52329 


+- 


1 


63555 


591 


26274 


+- 





27003 


44000000 


588 


91748 


+- 


1 


60674 


591 


20943 


+- 





26641 


45000000 


590 


68588 


+- 


1 


61247 


591 


19780 


+- 





26294 


46000000 


590 


26056 


+- 


1 


80667 


591 


17743 


+- 





26021 


47000000 


590 


89561 


+- 


1 


66101 


591 


17143 


+- 





25711 


48000000 


593 


22822 


+- 


1 


64885 


591 


21428 


+- 





25409 


49000000 


591 


98897 


+- 


1 


63064 


591 


23009 


+- 





25112 


50000000 


593 


47276 


+- 


1 


71196 


591 


27494 


+- 





24847 



Result : 



Number of weighted events = 50000000 
Average = 591.2749427045 fb 

Standard deviation = 0.2484667990 fb 
Maximal weight = 0.0287650976 fb 
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Subcontributions without naive QCD: 



Tree-level four-fermion cross section: 
Average = 627.6779997923 fb 

Standard deviation = 0.1075731739 fb 
Maximal weight = 0.0003787331 fb 

2->4 DPA part : 

Average = -6.1930315855 fb 

Standard deviation = 0.0030507290 fb 

Maximal weight = 0.0000240297 fb 

2->4 virtual+sof t+f inal-state collinear photon part: 

without LL (alpha) ISR 

Average = -67.4968488454 fb 

Standard deviation = 0.0117625527 fb 

Maximal weight = 0.0000250099 fb 

2->4 initial-state collinear photon part (e+) : 
without LL (alpha) ISR 
Average = -65.5094255283 fb 

Standard deviation = 0.0131360874 fb 
Maximal weight = 0.0000551979 fb 

2->4 initial-state collinear photon part (e-) : 
without LL (alpha) ISR 
Average = -65.4479754318 fb 

Standard deviation = 0.0131306505 fb 
Maximal weight = 0.0000523650 fb 

2->4 LL initial-state radiation up to 0(alpha~3): 
Average = -54.3690820517 fb 

Standard deviation = 0.0645784370 fb 
Maximal weight = 0.0003228727 fb 

2->5 Bremsstrahlung contribution (QED) : 
Average = 201.0338802635 fb 

Standard deviation = 0.2234350653 fb 
Maximal weight = 0.0277012770 fb 

2->5 Bremsstrahlung contribution (QCD) : 
Average = 0.0000000000 fb 

Standard deviation = 0.0000000000 fb 
Maximal weight = 0.0000000000 fb 
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7. events 


rejected 


in 


calculation 


of 


weight ( 


1) 


5 


3% 


7. events 


rejected 


in 


calculation 


of 


weight ( 


2) 


5 


27. 


7. events 


rejected 


in 


calculation 


of 


weight ( 


3) 


5 


a 


7. events 


rejected 


in 


calculation 


of 


weight ( 


4) 


6 


a 


7. events 


rejected 


in 


calculation 


of 


weight ( 


5) 


6 


2% 


7. events 


rejected 


in 


calculation 


of 


weight ( 


6) 


6 


07. 


% events 


rejected 


in 


calculation 


of 


weight ( 


7) 


34 


67. 


% events 


rejected 


in 


calculation 


of 


weight ( 


8) 





07. 



A. 2 Test run in the IB A mode for the processes e + e — > 4/7 

The second test run calculates the IB A for the process e + e~ — > ud/z~z/^7 in the subtraction 
branch for a CM energy of = 200 GeV. The ADLO cuts are applied, and the IBA W- 
boson width is calculated. On a DEC/ALPHA workstation the computing time is about 
2 hours. 

A. 2.1 Input file input iba4f a 



outputfile ! name of output file 

200d0 ! energy: CMF energy (in GeV) 

10000000 ! neventsw: number of weighted events 

2 ! smc: choice of MC branch: l(or 3): slicing 2 : subtraction 

! sborn4: include Born ee->4f : 0:no 1-3: yes 

1 ! sborn5: include Born ee->4f +photon: 0:no l:yes 

! sborng5: include Born ee->4f +gluon: 0:no l:yes 

1 ! sisr: include higher-order ISR: 0:no l:yes 

3 ! src: include radiative corrections: 0:no 1:DPA 2:IBA-4f 3:IBA-4fa 
1 ! scoultree: Coulomb singularity for ee->4f ,4f +ga: 0:no 1,2: yes 

3 ! qnf : Coulomb singularity for ee->4f (DPA) : 1,2, or 3 

1 ! qreal: neglect imaginary part of virt . corr. : 0:no l:yes 

2 ! qalp: choice of input -parameter scheme: 0,1, or 2 
2 ! qgw: calculate the W-boson width: 0:no 1-4: yes 

1 ! qprop: choice of width scheme: 0,1,2,3 or 4 

! ssigepem4: choice of diag. for Born ee->4f : 0:all 1-5: subsets 

! ssigepem5: choice of diag. for Born ee->4f+ga: 0:all 1-5: subsets 

! ssigepemg5: choice of diag. for Born ee->4f+gl: 0:all 1,5: subsets 

2 ! qqcd: include QCD radiative corr.: 0:no 1:CC03 2:naive 3:CC11 
! sqcdepem: include gluon-exch. diag. in Born: 0:no l:yes 2: only 
u ! fermion 3 

d ! anti-fermion 4 

mu ! fermion 5 

nu_mu ! anti-fermion 6 

OdO ! pp: degree of positron beam polarization [$-l$d0:ld0] 

OdO ! pm: degree of electron beam polarization [$-l$d0:ld0] 

! srecomb: recombination cuts: 0:no 1:TH 2: EXP 
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OdO 


! precomb(l): angular rec. cut between photon and beam 


OdO 


! precomb(2) : rec. cut on photon energy 


OdO 


! precomb(3): inv.-mass rec.(TH) or angular rec. cut for lept.(EXP) 


OdO 


' Drecombf4) * angular rec cut for auarks CEXP") 

• J- ^— ' W Will \J \ J. J * U-llp. U. -L- V_i_l_ -l_ W W • W 1-4. \J _J_ W J_ \J U.U.X Xl \ J — 1 ill / 


o 


1 Rrprnmbcr * crl n on TPromVii n at i on mt.s ■ O*no 1 * TH 9 * F.XP 

• O J_ v— ' v_> Will l_/ t-^ * EC. _L LX WXX J- \wf will _L XX CX O _L ull v_> IX \j O * W * XX W _L • XXX * XJiLl 


OdO 


! precombg(l) : rec. cut on gluon energy 


OdO 


! precombg (2) : inv . -mass (TH) or angular (EXP) recombination cut 


o 


l sat.pr ■ an nma Ions trinle can ere conn! i ncrs (TClC,} * * n o 1 *vpr 

• kJ <—L L/Clj * qXlVllig.X. W UkJ OX. -L. _l_ V-^ cLLLUcL v— ' WW 1-4. !_/ -1 LllcL >_> \ 1 Uu/ • \J t Uu J- • y \^ kJ 


OdO 


! TGC Delta g_l~A 


OdO 


! TGC Delta g_l~Z 


OdO 


! TGC Delta kappa~A 


OdO 


1 TGC Delta kaDDa~Z 


OdO 


! TGC lambda~A 


OdO 


! TGC lambda" Z 


OdO 


! TGC g_4"A 


OdO 


! TGC g 4~Z 


OdO 


' TGC e- 5~A 

• X W \J EL w 


OdO 


' TGC e- 5~Z 

. X W \_* EL W i_j 


OdO 


• TGC tilde kar>r>a~A 

. x w \_* <-<_l,_l.w.w Xi.cj.www. n 


OdO 


• TGC tilde kar>na~Z 

X \J \J U_L_LW.W HUk/L/U 


OdO 


! TGC tilde lambda~A 


OdO 


I TGC tilde lambda~Z 


OdO 


! TGC f 4~A 


OdO 


! TGC f 4~Z 


OdO 


! TGC f_5~A 


OdO 


! TGC f_5~Z 


OdO 


! TGC h_l~A 


OdO 


! TGC h_l~Z 


OdO 


! TGC h_3~A 


OdO 


! TGC h_3~Z 





! qaqgc : anomalous quartic gauge couplings (QGC) : 0:no l:yes 


OdO 


! QGC a_0/Lambda~2 


OdO 


! QGC a_c/Lambda"2 


OdO 


! QGC a_n/Lambda"2 


OdO 


! QGC tilde a_0/Lambda~2 


OdO 


! QGC tilde a_n/Lambda~2 


1 


! scuts: separation cuts: 0:no l,2:default(ADL0,LC) 10,ll:input 



A. 2.2 Output file 



smc= 2: Subtraction-method branch of RacoonWW 



Input parameters: 



CMF energy = 200.00000 GeV, Number of events = 10000000, 
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alpha (0) 
GF 
MW 
GW 
me 
mu 
md 



= 1/ 137.0359895, alpha(MZ) = 1/128.88700, alpha_s = 0.11900, 

MH = 150.00000, 



. 1166370E-04, 
80.35000, 

2.09436, 
.51099907E-03, 

0.00485, 

0.00485, 



MZ = 91.18670, 

GZ = 2.49471, 

mmu = 0.105658389, mtau = 1.77705, 

mc = 1.55000 mt = 174.17000, 

ms = 0.15000, mb = 4.50000. 



Effective branching ratios: 

leptonic BR = 0.32512, hadronic BR = 0.67488, total BR = 1.00000 
Process: anti-e e -> u anti-d mu anti-nu_mu + photon 



pp= 
pm= 
qalp= 
qgw= 
qprop= 

sborn4= 
src= 
scoultree= 
ssigepem5= 
qqcd= 
sisr= 



scuts= 1: 

energy (3) > 

energy (4) > 

energy (5) > 

energy (7) > 

angle (1,5) > 

angle (1,7) > 

angle (2, 5) > 

angle (2, 7) > 

angle (3, 5) > 

angle (3, 7) > 

angle (4, 5) > 

angle (4, 7) > 

angle (5, 7) > 

mass (3, 4) > 

events 
1000000 
2000000 
3000000 
4000000 



0.0: degree of positron beam polarization. 
0.0: degree of electron beam polarization. 
2: GF-parametrization scheme. 
2: IBA W-boson width calculated. 
1: constant width. 

0: tree-level process ee -> 4f NOT included. 
3: improved Born approximation ee->4f +photon. 
1: Coulomb sing, to ee->4f +photon included. 
0: all electroweak diagrams included. 
2: naive QCD corrections included. 

1: initial-state radiation up to order alpha~3 included. 

with separation cuts: 
3.00000 GeV 
3.00000 GeV 
1.00000 GeV 
0.10000 GeV 
10.00000 deg 
1.00000 deg 
10.00000 deg 



,00000 deg 
,00000 deg 
,00000 deg 
.00000 deg 
.00000 deg 
.00000 deg 
.00000 GeV 



intermediary results 
185.52994 +- 1.33873 
184.94817 +- 1.21671 
187.67723 +- 1.22280 
188.12389 +- 1.29820 



preliminary results 
185.52994 +- 1.33873 
185.23906 +- 0.90451 
186.05178 +- 0.72784 
186.56981 +- 0.63508 
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5000000 
6000000 
7000000 
8000000 
9000000 
10000000 



186.24278 +- 
184.67305 +- 
186.27700 +- 
188.42698 +- 
188.10039 +- 
187.29950 +- 



1.21187 
1.27906 
1.29118 
1.29575 
1.28908 
1.27811 



186.50440 +- 
186.19918 +- 
186.21029 +- 
186.48738 +- 
186.66660 +- 
186.72989 +- 



0.56291 
0.51526 
0.47862 
0.44903 
0.42406 
0.40248 



Result : 



Number of weighted events = 10000000 

Average = 186.7298922322 fb 

Standard deviation = 0.4024824609 fb 

Maximal weight = 0.0452822218 fb 



Subcontributions without naive QCD: 



Tree-level f our-f ermion cross section: 
Average = 0.0000000000 fb 

Standard deviation = 0.0000000000 fb 
Maximal weight = 0.0000000000 fb 
(event = 0, channel = 0) 



2->4 part of radiative corrections without (alpha) LL ISR: 
Average = 0.0000000000 fb 

Standard deviation = 0.0000000000 fb 
Maximal weight = 0.0000000000 fb 
(event = 0, channel = 0) 



Initial-state radiation up to 0(alpha~3): 
Average = 0.0000000000 fb 

Standard deviation = 0.0000000000 fb 
Maximal weight = 0.0000000000 fb 
(event = 0, channel = 0) 



Improved Born Approximation ee->4f +photon: 
Average = 179.9149173069 fb 

Standard deviation = 0.3877932869 fb 
Maximal weight = 0.0436295823 fb 
(event = 5915698, channel = 50) 



% events 


rejected 


in 


calculation 


of 


weight ( 


1) 





0% 


% events 


rejected 


in 


calculation 


of 


weight ( 


2) 





0% 


% events 


rejected 


in 


calculation 


of 


weight ( 


3) 





0% 


% events 


rejected 


in 


calculation 


of 


weight ( 


4) 





0% 


% events 


rejected 


in 


calculation 


of 


weight ( 


5) 





0% 


% events 


rejected 


in 


calculation 


of 


weight ( 


6) 


10 


2% 


% events 


rejected 


in 


calculation 


of 


weight ( 


7) 





0% 
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°/ events rejected in calculation of weight(29): 0.0% 
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